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Abstract 

This  research  investigates  a  microstrip  interdigitated  electrode  structure  used  in  integrated 
circuit  sensors  and  presents  an  electrical  model  for  the  device.  The  structure  is  assumed  to:  (1) 
have  a  periodic  arrangement  of  interdigitated  electrodes,  (2)  have  electrodes  of  sufTuient  length  to 
ignore  effects  created  by  the  ends  of  the  structure,  (3)  be  operating  at  a  sufficiently  low  frequency, 
.vhere  the  potential  is  essentially  constant  along  the  electrode’s  length  and  where  the  substrate  is 
essentially  a  conductor,  and  finally  (4)  have  five  layers  of  homogeneous,  isotropic,  non-magnetic 
material  in  the  plane  transverse  to  the  structure’s  length.  Within  these  constraints,  the  problem 
focuses  on  solving  for  the  potential  distribution  using  Laplace’s  equation. 

The  electrical  model  developed  consists  of  three  elements — a  capacitance  from  each  of  two 
electrodes  to  ground,  a  capacitance  between  the  electrodes,  and  a  resistance  through  the  sensor 
coating.  The  values  for  the  individual  elements  are  found  from  the  potential  solution.  Potential 
functions,  in  the  form  of  infinite  Fourier  series,  are  developed  for  each  layer  witliin  a  periodic 
cell.  The  coefficients  of  the  Fourier  series  are  determined  by  matching  boundary  conditions  across 
interfaces  between  the  material  layers.  These  coefficients  are  estimated  numerically  by  truncating 
the  infinite  series,  and  solving  a  linear  system  of  equations. 

Results  are  presented  for  a  variety  of  cell  parameters  (cell  dimensions  and  material  properties), 
which  illustrate  the  dependence  of  the  individual  elements  on  these  cell  parameters.  Specifically, 
the  results  indicate  that  all  model  elements  are  strongly  influenced  by  the  electrical  properties  of  the 
sensor  coating;  and  therefore,  all  elements  play  a  major  role  in  the  sensing  action  of  the  device.  The 
frequency  dependence  of  the  model  parameters  are  also  developed,  using  a  complex  permittivity. 
Results  are  presented  for  various  material  parameters,  which  illustrate  the  high  and  low  frequency 
limits  of  all  model  elements  with  a  transition  centered  at  the  frequency  where  the  sensor  coating 
changes  from  a  conductor  to  a  lossy  dielectric. 
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RIGOROUS  CALCULATIONS  OF  LUMPED  ELECTRICAL 
PARAMETERS  FOR  A  MICROSTRIP  INTERDIGITATED 
ELECTRODE  STRUCTURE 


I.  Introduction 


1.1  Purpose 

This  thesis  investigates  the  microstrip  interdigitated  electrode  structure  used  in  integrated 
circuit  sensors  and  presents  an  electrical  model  for  the  device.  The  electrical  parameters  of  interest 
are  developed  using  the  potential  distribution  within  the  device. 

1.2  Organization 

The  remainder  of  this  chapter  contains  definitions  of  key  terms  and  summaries  of  previ¬ 
ous  research  which  motivated  this  thesis  topic.  Chapter  II  presents  the  specific  thesis  problem 
and  reviews  related  research  investigating  similar  topics.  Chapter  III  provides  an  analysis  of  the 
interdigitated  electrode  sensor  structure,  which  leads  to  a  simplified  mathematical  model.  This 
analysis  also  introduces  an  approximation  under  which  the  potential  distribution  satisfies  Laplace’s 
equation  in  the  transverse  plane;  it  also  provides  the  mathematical  conditions  that  this  transverse 
pntpntial  distribution  must  satisfy.  Chapter  IV  presents  the  solution  for  the  capacitive  elements, 
while  Chapter  V  presents  the  solution  for  the  resisitive  element.  Numerical  results  of  the  research 
are  presented  for  a  particular  device  in  Chapter  V’l,  including  graphs  of  the  model  elements  with 
respect  to  variations  of  the  physical  and  material  aspects  of  the  structure.  In  addition,  constant 
voltage  contour  plots  of  the  structure  arc  presented  to  illustrate  the  solution.  Chapter  VII  presents 
a  frequency  analysis  of  the  model  using  complex  permittivity,  which  displays  the  effects  of  the 
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sensor  coating’s  permittivity  and  conductivity  on  the,  now,  frequency-dependent  lumped  elements. 
The  thesis  concludes  with  a  brief  summary  and  recommendations  for  further  research. 

1.3  Definiiions 

The  key  terms  related  to  this  research  effort  are  microstrip  transmission  line,  interdigitated 
electrode  structure,  and  interdigitated  electrode  sensor  structure. 

A  microstrip  is  a  conducting  metal  strip  used  in  electronic  circuits  to  convey  electrical  energy 
from  one  circuit  component  to  another.  A  layer  of  insulating  dielectric  material  separates  the  strip 
of  metal  from  a  substrate  having  bulk  dielectric  and  conductive  properties  (see  Figure  1.1).  A 
ground  plane  is  located  below  the  substrate. 


Figure  1.1.  Microstrip  transmission  line 


.4n  interdigitated  electrode  structure  is  a  configuration  similar  to  two  forks  with  their  tines 
interlaced.  The  microstrip  interdigitated  electrode  structure  is  a  configuration  of  interlaced  mi- 


crostrips  (fingers),  where  each  microstrip  segment  is  located  between  two  microstrips  from  the 
other  electrode  (see  Figure  1.2). 


Figure  1.2.  Microstrip  interdigitated  electrode  structure 

An  interdigitated  electrode  sensor  structure  uses  the  microstrip  interdigitated  electrode  struc¬ 
ture  with  a  special  coating  applied  over  the  structure  (see  Figure  1.3).  The  coating  material  has 
variable  electrical  properties,  which  are  dependent  on  external  events  and  which  change  the  elec¬ 
trical  performance  of  the  structure.  The  sensor’s  fundamental  purpose  is  to  detect  these  changes. 

1.4  Background 

Integrated  circuit  sensors  are  miniaturized  versions  of  the  interdigitated  electrode  sensor  struc¬ 
ture  along  with  support  electronics  necessary  to  monitor  its  electrical  performance.  Two  such  de¬ 
vices  are  offered  as  examples — a  chemically-sensitive  field-effect  transistor  (CHEMFET)  and  a  resin 
cure  sensor. 
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Figure  1.3.  Microstrip  interdigitated  electrode  sensor  structure 


1.41  CHEMFET.  Wiseman  investigated  two  different  (small  and  large)  interdigitated  elec¬ 
trode  sensor  structures  utilizing  the  CHEMFET  configuration  [18].  He  coated  the  interdigitated 
structure  with  a  special  polymer  (copper  phthalocyanine)  whose  electrical  properties  are  sensitive  to 
certain  environmentally-sensitive  gas  contaminants  (nitrogen  dioxide  and  diisopropyl  methylphos- 
phonate).  The  large  interdigitated  structure  had  29  inner  fingers  and  30  outer  fingers.  These  fingers 
were  7.5  ftm  wide,  3.792  mm  long,  and  they  were  separated  from  each  other  by  9.0  ^m.  The  small 
interdigitated  structure  had  15  inner  fingers  and  16  outer  fingers.  These  fingers  were  6.0  /im  wide, 
3.792  mm,  and  were  similarly  separated  from  each  other  by  9.0  /im. 

Wiseman  did  not  mathematically  analyze  his  interdigitated  structure  (from  an  electromag¬ 
netic  perspective),  and  the  dimensions  were  chosen  based  upon  fabrication  design  rules.  Wiseman 
cited  work  by  Lee  (section  1.4.2)  as  validation  for  using  the  interdigitated  structure.  Wiseman  rec- 
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ommended  further  research  regarding  the  interdigitated  structure,  with  the  goal  of  mathematically 
predicting  a  physical  configuration  which  offers  the  greatest  electrical  response  to  a  given  variation 
in  the  coating  material’s  electrical  properties,  thereby  creating  an  optimum  sensing  capability. 

1.4-2  Epoxy  Resin  Cure  Sensor.  Lee  presented  an  epoxy  resin  cure  sensor  using  an  inter¬ 
digitated  structure  [14].  He  coated  the  interdigitated  structure  with  epoxy  resin  whose  electrical 
properties  vary  depending  on  the  epoxy’s  cure  state.  Lee’s  interdigitated  structure  had  nine  inner 
and  ten  outer  fingers.  The  fingers  were  12.5  /jm  wide,  5  mm  long,  and  were  separated  from  each 
other  by  12.5  ^m. 

Lee  performed  a  detailed  analysis  of  his  interdigitated  structure.  He  developed  a  circuit  model 
of  the  structure  using  a  finite  difference  numerical  solution  to  determine  the  potential  distribution. 
Using  the  finite  difference  method,  unknown  voltages  were  determined  at  discrete  points  throughout 
the  structure.  This  led  to  a  large  number  of  unknowns  (10,000  for  a  100  by  100  point  grid). 
Since  the  potential  distribution  is  valid  only  for  a  specific  physical  configuration,  these  unknowns 
were  determined  for  each  structural  arrangement  under  investigation.  While  the  finite  difference 
technique  yields  valid  results,  it  can  be  advantageous  to  solve  for  the  potential  distribution  with 
a  method  requiring  fewer  unknowns.  In  addition,  Lee’s  analysis  assumed  the  microstrips  were 
infinitely  thin,  which  limited  his  results.  For  future  research,  Lee  recommended  a  full  analysis  of 
the  interdigitated  structure  using  microstrips  with  finite  height. 

1.4.3  Background  Analysis.  In  each  of  the  integrated  sensor  circuits  described,  most  of  the 
interdigitated  electrode  structure  dimensions  were  influenced  by  fabrication  limitations.  While 
each  of  these  research  efforts  demonstrated  the  effectiveness  of  the  interdigitated  electrode  sensor 
structure,  each  recommended  further  research  to  investigate  a  physical  configuration  which  would 
yield  the  best  sensor,  or  simply  have  the  greatest  electrical  response  to  any  changes  in  the  applied 
coating’s  material  properties. 
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II.  Problem  Statement 


2.1  Research  Objective:  Electrical  Model 

To  adequately  predict  an  optimum  design  for  the  interdigitated  electrode  sensor  structure, 
a  lumped  element  electrical  model  of  the  device  is  proposed.  The  problem  solved  in  this  thesis  is 
the  development  of  lumped  electrical  parameters  to  describe  the  device.  This  chapter  describes  the 
specific  electrical  parameters  being  addressed  and  reviews  the  literature  concerning  the  methods 
utilized  to  determine  similar  parameters. 

2.2  Identification  of  Lumped  Electrical  Parameters 

The  individual  parameters  under  investigation  are  described  below.  They  are  a  capacitance 
associated  with  the  substrate  and  insulator,  resistance  of  the  substrate,  capacitance  between  elec¬ 
trodes,  and  resistance  between  electrodes  (Figure  2.1). 

Capacitance  Associated  with  the  Substrate  and  Insulator.  For  an  applied  voltage, 
there  is  a  corresponding  charge  distribution  on  each  finger.  This  charge  can  be  used  to  define 
a  capacitance  value  between  each  finger  and  the  ground  plane  below.  The  total  capacitance  is 
the  combination  of  a  capacitance  associated  with  the  substrate  (labeled  C,  in  Figure  2.1)  and  a 
capacitance  associated  with  the  insulator  (labeled  C,-  in  Figure  2.1). 

Resistance  of  the  Substrate.  The  substrate  is  composed  of  a  material  having  bulk 
conductive  properties.  Thus,  the  substrate  will  provide  some  resistance  between  the  oxide  and  the 
ground  plane.  This  resistance  is  labeled  R,  in  Figure  2.1. 

Capacitance  Between  Electrodes.  A  capacitance  between  the  electrodes  can  be 
defined  based  upon  the  charge  and  potential  differences  between  the  two  electrodes.  This  capaci- 
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Figure  2.1.  Microstrip  interdigitated  electrode  sensor  structure  parameters 


tance  value  varies  due  to  the  electrical  properties  of  the  sensor  coating.  The  capacitance  between 
electrodes  is  labeled  Ce  in  Figure  2.1. 

Resistance  Between  the  Electrodes.  Because  the  sensor  coating  is  partially  con¬ 
ductive,  there  is  a  finite  resistance  between  the  two  electrodes.  This  resistance  typically  varies  with 
respect  to  the  chemical  properties  and  state  of  the  sensor  coating,  and  it  is  labeled  Rp  in  Figure  2.1. 

2.3  Literature  Review  of  Microstrip  Parameters 

Several  authors  have  investigated  microstrip  parameters  in  various  configurations.  This  sec¬ 
tion  reviews  microstrip  interdigitated  electrode  structures  and  miscellaneous  microstrip  configura¬ 
tions. 
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2.3.1  Microstnp  'nterdigitated  Electrode  Structures.  The  literature  concerning  other  inter- 
digitated  electrode  structures  typically  offers  analysis  of  the  fingers  in  a  specific  configuration. 
Unfortunately,  the  interdigitated  electrode  sensor  structure  is  unique  because  of  the  partially  con¬ 
ductive  coating  covering  the  electrode  fingers.  However,  the  approach  and  analysis  of  similar 
structures  without  sensor  coatings  provides  insight  into  ihe  techniques  that  can  be  employed  to 
analyze  the  interdigitated  electrode  sensor  structure. 

Interdigital  Capacitor.  Alley  analyzed  a  thin  film  interdigital  capacitor  [2].  His 
interdigitated  structure  contained  an  even  number  of  fingers  (equal  number  per  electrode)  and  the 
fingers  had  equal  width  and  separation.  Alley  detemined  the  capacitance  of  his  structure  using  even 
and  odd  mode  capacitance  formulas  for  coupled  lines  developed  by  Smith  [15].  While  Alley  did 
consider  possible  resistance  between  the  electrodes  (due  to  a  lossy  dielectric  beneath  the  fingers), 
his  analysis  was  limited  to  infinitely  thin  electrodes  (no  finite  height). 

Interdigitated  Microstrip  Parameters.  Dalby  presented  interdigitated  microstrip 
circuit  parameters  for  both  single  and  coupled  microstrips  [5].  His  physical  arrangement  included 
an  arbitrary  material  layer  above  the  microstrips  and  arbitrary  dimensions  for  all  finger  features. 
Dalby  presented  parameter  formulas  for  both  thin  and  thick  single  microstrips,  ^ls  well  as  odd  and 
even  mode  parameter  formulas  for  the  thin  coupled  microstrips.  Unfortunately,  these  formulas  were 
presented  without  derivation  or  reference,  which  limited  this  reader’s  ability  to  modify  his  results 
for  other  configurations  or  verify  the  accuracy  of  his  formulas. 

Smith  developed  formula  for  even  and  odd  mode  capacitances  for  coupled  microstrips  using 
conformal  mapping  techniques  (15).  As  with  Dalby’s  work.  Smith  considered  arbitrary  material 
layers  above  and  below  the  microstrips  with  arbitrary  dimensions  for  the  coupled  strips.  While 
Smith  rigorously  presented  the  upper  and  lower  bounds  to  the  capacitance  values,  his  results  do 
not  apply  to  this  specific  problem  because  he  employed  a  shielding  enclosure  and  microstrips  of 
infinitesimal  height. 
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2.3.2  Miscellaneous  Microstrxp  Configurations.  Aside  from  interdigiiated  structures,  many 
articles  have  been  published  concerning  microstrip  transmission  lines  that  may  be  useful  in  analyz¬ 
ing  the  interdigitated  electrode  sensor  structure.  In  all  of  the  articles  discussed  below,  the  electrical 
parameters  were  calculated  by  first  determining  the  electric  field  associated  with  a  specified  physical 
configuration.  While  none  of  the  arrangements  match  the  interdigitated  electrode  sensor  structure 
exactly,  the  generalized  method  provides  a  logical  approach  for  developing  a  solution  to  this  thesis 
problem  (seek  the  electric  field  associated  with  the  interdigitated  electrode  sensor  structure). 

Microstrips  in  Various  Configurations.  Yamashita  and  Atsuki  analyzed  shielded 
stripline  within  three  layers  of  material  using  a  Green’s  function  and  a  variational  method  [19]. 
They  developed  expressions  for  the  line  capacitance  and  other  strip  parameters  (characteristic 
impedance,  wavelength,  and  wave  velocity).  While  they  considered  the  thickness  of  the  stripline 
conductors,  they  did  not  address  coupled  strips.  Because  of  this  limitation  and  their  shielded 
enclosure,  their  results  did  not  apply  to  the  interdigitated  electrode  sensor  structure. 

Kitazawa  and  Hayashi  analyzed  two-line  [10]  and  three-line  [11]  coupled,  thin  striplines  em¬ 
bedded  in  three  layers  of  anisotropic  material  with  ground  planes  above  and  below  these  layers. 
They  developed  a  capacitance  matrix  using  a  Green’s  function  and  a  variational  method.  Ex¬ 
pressions  for  the  line  capacitance  and  other  line  parameters  (charateristic  admittance,  impedance 
matrix,  and  inductance)  were  presented.  Since  they  neglected  strip  height  and  used  an  upper 
ground  plane,  their  results  are  not  adaptable  to  this  thesis  problem. 

Krage  and  Haddad  investigated  shielded  single  and  coupled  thin  microstrips  embedded  in 
two  and  three  layers  of  material  [12].  They  solved  for  the  electric  and  magnetic  fields,  using 
a  Fourier  series  representation  for  the  fields,  and  determined  the  Fourier  coefficients  from  the 
boundary  conditions  between  layers.  They  presented  expressions  for  the  line  capacitance  and  other 
parameters  (characteristic  impedance  and  line  inductance)  bEised  upon  the  field  solutions.  Since 
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the  microstrips  were  shielded  and  strip  height  was  neglected,  their  results  do  not  apply  to  the 
interdigitated  electrode  sensor  structure. 

Judd  and  others  analyzed  shielded  single  and  coupled  thick  microstrips  embedded  in  two  layers 
of  material  [8].  A  Fourier  series  representation  of  the  potential  in  each  layer  was  presented,  and  the 
Fourier  coefficients  were  determined  using  the  boundary  conditions  between  layers.  They  presented 
an  expression  for  the  capacitance  based  upon  the  potential  solution,  and  the  capacitance  was 
estimated  by  truncating  the  infinite  series.  These  specific  results  did  not  apply  to  the  thesis  problem, 
because  a  shielded  enclosure  was  incorporated  and  only  two  layers  of  material  was  considered. 
However,  this  method  did  seem  suitable  for  analyzing  the  interdigitated  electrode  sensor  structure. 

Hill  and  others  investigated  coupled  microstrips  and  triplates  in  a  single  medium  [6].  They 
solved  for  the  electric  field  within  the  structure  using  an  infinite  series  form  of  a  Green’s  function. 
An  expression  for  the  capacitance  was  then  developed  from  the  field  solution.  The  specific  results 
do  not  apply  to  the  interdigitated  electrode  sensor  structure,  because  their  analysis  incorporated 
only  one  material  layer. 

General  Multi-Conductor  Capacitance  Matrix.  Weeks  presented  a  generalized 
method  for  determining  a  capacitance  matrix  for  multi-conductor  transmission  lines  [17].  With 
this  technique,  the  conductors  are  set  to  specific  potential  levels,  which  isolates  each  capacitance 
term  and  allows  each  capacitive  element  to  be  determined  separately.  Kitazawa  used  this  method 
for  the  Ccise  of  a  multi-microstrip  configuration  [9].  While  these  microstrips  were  not  arranged  in 
an  interdigitated  configuration,  Kitazawa  proved  that  Weeks’  technique  is  also  suitable  for  ana¬ 
lyzing  microstrips.  Kitaizawa  developed  expressions  for  the  electric  and  magnetic  fields  within  his 
arbitrary  arrangement  of  equal  height  microstrips  using  a  dyadic  Green’s  function.  Once  the  fields 
were  known,  Kitazawa  used  Weeks’  method  to  determine  all  resulting  capacitance  parameters.  Be¬ 
cause  Kitazawa  did  not  consider  conductive  materials,  his  specific  results  do  not  apply  to  the  thesis 
problem. 
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2.4  Model  Elements  from  the  Electric  Field 

The  importance  of  the  electric  field  was  illustrated  throughout  the  previous  articles.  Once 
the  electric  field  was  determined  for  a  particular  microstrip  structure,  the  electrical  parameters  of 
interest  were  found  using  fundamental  relationships,  such  as  Gauss’  law  and  Ohm’s  law. 


Capacitive  Elements.  The  charge,  Q,  on  the  microstrip  electrode  is  found  using 
Gauss’  law  by  integrating  the  electric  flux  density  over  a  suitable  surface,  7,  enclosing  the  electrode. 
Capacitance,  C,  is  then  found  from  the  relationship  Q  =  CV ,  where  V  is  the  voltage  potential. 
Using  a  forcing  potential  of  1  volt  yields 


(2.1) 


where  e  is  the  permittivity  of  the  material,  is  the  electric  field  normal  to  the  surface,  and  dl  is 
an  element  of  surface  area. 


Resistive  Elements.  Within  conductive  regions,  the  current  density,  J,  is  found 
using  the  relationship  J  =  <tE,  where  cr  is  the  conductivity  of  the  material  and  E  is  the  electric 
field.  The  current  through  the  conductive  region  is  found  by  integrating  the  current  density  over 
a  suitable  surface  (7),  through  all  of  the  current  paths.  The  resistance  is  then  found  using  Ohm’s 
law,  V  =  IR,  where  I  is  the  current  and  R  is  the  resistance,  or  equivalently  G  =  I/V .  where  G  is 
the  conductance  (G  =  l/R).  Using  a  forcing  potential  of  1  volt  yields 

G  =  I  J„dl  =  f  <TEr,dl  (2.2) 

Jy  Jy 

where  and  £„  are  the  current  density  and  electric  field,  respectively,  normal  to  the  surface. 

Therefore,  to  develop  the  electrical  model  elements,  we  seek  a  solution  to  the  electric  field 
within  the  microstrip  interdigitated  electrode  sensor  structure. 
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in.  Problem  Analysis 


It  was  shown  in  the  previous  chapter  that  the  electrical  parameters  of  interest  can  be  deter¬ 
mined  from  the  electric  field  within  the  interdigitated  electrode  structure.  Since  the  electric  field 
can  be  found  directly  from  the  potential  distribution,  the  thesis  concentrates  on  solving  for  the 
potential  distribution  within  the  interdigitated  electrode  sensor  structure.  This  research  further 
concentrates  on  the  plane  transverse  to  the  length  of  the  structure,  since  the  length  is  the  largest 
single  dimension  (an  electrode  finger  is  much  longer  than  its  width  or  height).  Therefore  the  effects 
at  the  structure’s  endpoints  (the  bus  where  the  appropriate  digits  are  joined)  are  neglected.  Ac¬ 
cording  to  Alley  [2:1029],  neglecting  end-effects  introduces  only  a  small  error  (less  than  2%)  in  the 
capacitance  values.  F.ncl-effects  are  also  ignored  in  this  thesis  for  the  resistive  element,  although  no 
analysis  is  provided,  nor  is  any  reference,  to  support  this  approximation.  Towards  the  objective  of 
solving  for  the  potential,  the  physical  structure  must  be  analyzed  to  determine  the  mathematical 
conditions  that  the  potential  must  satisfy. 

In  section  3,1.  the  periodic  nature  of  the  interdigitated  electrode  sensor  structure  is  addressed 
and  a  corresponding  periodic  cell  is  created  to  represent  the  interdigitated  structure,  while  the 
material  layers  within  the  cell  are  defined  in  section  3.2.  In  section  3.3,  a  relationship  is  developed 
between  the  electric  field  and  the  transverse  potential,  using  a  transverse  electromagnetic  (TEM) 
wave  analysis.  An  upper  frequency  bound  for  the  analysis  i.s  then  established  in  section  3.4,  which 
reduces  the  electric  field  to  simply  the  gradient  of  the  transverse  potential.  In  section  3.5,  this 
potential  is  shown  to  satisfy  Laplace's  equation.  Section  3  6  discusses  the  boundary  conditions 
between  the  layers,  which  lead  to  specific  boundary  conditions  that  a  high  and  low  frequency 
potential  solution  must  satisfy.  Finally,  in  section  3.7,  a  summary  is  presented  to  consolidate  the 
analysis. 
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3.]  Periodic  Structure 


As  viewed  from  the  transverse  plane  (transverse  to  the  length  of  the  electrodes),  the  interdig- 
itated  electrode  sensor  structure  possesses  a  repetitive  pattern  consisting  of  two  adjacent  electrcde 
fingers  (refer  to  Figure  1.3).  Since  the  interdigitatcd  electrode  sensor  structure  contains  many 
fingers  (59  in  the  CIIEMFET),  this  pattern  repeatr  itself  30  times.  Hence,  a  periodic  cell  can 
be  created  to  represent  the  interdigitatcd  electrode  structure.  This  approximation  introduces  an 
error  for  the  outer  few  sets  of  fingers,  but  fortunately,  this  perturbation  only  represents  a  small 
portion  of  the  total  structure.  The  arrangement  consisting  of  one-half  of  each  finger  is  considered 
the  periodic  cell  (see  Figure  3.1).  This  periodic  condition  requires  a  vanishing  normal  component 
ot  the  electric  field  on  the  vertical  walls  of  the  cell.  Considering  the  vertical  direction  to  be  the 
y-axis,  this  situation  implies  that  at  the  cell  walls 


Ex  =  0 

where  Ex  is  the  x  component  of  the  electric  field,  or  alternately, 

where  V  is  the  voltage,  or  potential. 


(3.1) 


(3.2) 


3.2  Five  Layers  of  Material 

There  are  five  layers  of  material  present  in  the  periodic  cell.  The  first  layer  is  the  substrate 
material,  wh  ch  is  directly  above  the  ground  plane.  The  second  layer  is  the  insulator  material,  which 
is  located  beneath  the  interdigitated  electrodes.  Tlie  third  layer  represents  the  region  between  the 
fingers.  Due  to  the  integrated  circuit  fabrication  process,  this  layer  is  part  insulti.;  (or  passivation) 
material,  part  sensor  coating,  and  part  air.  For  simplification,  this  layer  will  be  designated  insulator, 
but  the  model  will  remain  flexible  and  allow  arbitrary  material  parameters  in  a  homogeneous  layer 
3.  The  fourth  layer  is  the  sensor  coating,  while  the  fifth  layer  is  air.  Each  region  has  an  associated 
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Figure  3.1.  Periodic  cell  of  interdigitated  electrode  sensor  structure 


permittivity  (f),  while  the  substrate  and  sensor  coating  possess  finite  conductivity  (a)  as  well. 
The  permeability  of  each  region  is  considered  equal  to  the  permeability  of  free  space  These 

layers  are  also  identified  in  Figure  3.1.  To  simplify  the  model,  each  material  layer  is  considered 
homogeneous  and  isotropic;  that  is,  the  material  and  electrical  properties  are  uniform  throughout 
their  bulk  material. 

3.3  Transverse  Electomagneiic  Wave  in  the  Transverse  Plane 

In  solving  for  the  electric  field  within  the  interdigitated  electrode  sensor  structure,  it  is  nec¬ 
essary  to  establish  a  relationship  between  the  structure’s  electric  field  and  the  potential  within  the 
periodic  cell.  Towards  this  objective,  the  electrodes  are  considered  perfect  conductors,  and  all  the 
current  within  the  structure  is  considered  to  be  along  the  length  of  the  structure,  thereby  assuming 
transverse  electromagnetic  (TEM)  fields  within  the  microstrip  structure.  While  the  metal  may  not 
be  ideal  and  some  transverse  current  can  be  expected,  these  assumptions  provide  a  starting  point 
for  the  analysis  and  provide  a  basis  to  establish  a  relationship  bet\/een  the  potential  and  the  elec¬ 
tric  field.  The  following  analysis  summarizes  the  discussion  in  Collins  [3:67-70].  Let  z  denote  the 
direction  which  aligns  with  the  length  of  the  interdigitated  structure,  with  y  denoting  the  vertical 
direction  (with  respect  to  the  periodic  cell),  and  x  denoting  the  horizonal  direction  (with  respect 
to  the  periodic  cell).  The  basic  TEM  premise  is 

E,  =  H,  =  0  (3.3) 

where  Ei  and  Hz  are  the  z  components  of  the  electric  and  magnetic  fields,  respectively.  Equation 
3.3  implies  a  perfect  conductor,  and  all  the  current  is  in  the  z  direction.  The  corresponding  time- 
harmonic  Maxwell  equations,  for  exp(jw<)  time  dependence,  are 


V  X  E, 

(3.4) 

V  X  H, 

=  ju/eEt 

(3.5) 
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where  Ej  and  H,  are  the  transverse  electric  and  magnetic  fields,  respectively,  uj  is  the  angular 
(excitation)  frequency,  n  and  (  are  the  material’s  permeability  and  permittivity,  respectively,  and 
j  =  \/^.  Since 


V  =  X 


d 

dx 


d 

+  y-7r  + 


dz 


(3.6) 


where  i,  y,  and  i  are  the  unit  vectors  in  the  x,  y,  and  z  directions,  respectively.  Let 


V, 

V. 


d  d 


Then,  we  may  write 


V  =  V,  +  V, 


(3.7) 

(3.8) 


(3.9) 


Accordingly,  Maxwell’s  equations  can  be  expressed  as 


(V,  X  E.)  +  (V,  X  E,)  =  -juyn,  (3.10) 

(V,  X  H,)  +  (V,  X  H.)  =  jwcE,  (3.11) 


Using  Equation  3.3,  the  curl  of  the  transverse  electric  and  magnetic  fields  are 


V  X  E,  = 


dEy  dEi 


dx  dy 


(3.121 


and 


V  X  H, 


dHy  _  ^ 

dx  dy 


(3.13) 


Thus 


V(  X  E(  =  0 


(3.14) 
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V,  X  H, 


0 


(3.15) 


V;  X  Et  = 

(3.16) 

X  H/  =  jiuitEj 

(3.17) 

Since  Vt  x  E(  =  0,  Equations  3.14,  3.15,  3.16,  and  3.17  will  be  satisified  by 

E,  =  3i(z)V,0(r,y) 

(3.18) 

H«  =  y2(-)V,V'(i,y) 

(3.19) 

where  yi,  g^,  <t>,  and  t/’  are  scalar  functions  which  need  to  be  determined.  From 

Equation  3.16 

H,  =  X  E, 

JW/i 

(3.20) 

Using  Equation  3.20  in  3.17  yields 

V;  X  X  E,  =  w^/jcEf 

(3.21) 

Using  the  identity 

V,  X  V,  xE,  =  V,(V,  .E,)-V2e, 

(3.22) 

and  the  fact  that 

V,  ■  E,  =  0 

(3.23) 

yields 

X  X  E,  =  -V^E, 

(3.24) 

Using  Equation  3.24  in  3.21  yields 

(3.25) 
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where  and  k  is  the  wavenumber  associated  with  the  material.  Using  Equation  3.18  in 

3.25  yields 

y)  +  k-g,{z)V,4>(^,  y)  =  0  (3.26) 

where  the  prime  symbol  denotes  differentiation  with  respect  to  z.  For  a  non-trivial  transverse  field, 
this  implies 

=  0  (3.27) 

Hence 

ffi(^)  =  A.±  €xp(±jkz)  (3.28) 

where  A±  is  an  arbitrary  constant.  Using  Equation  3.28  in  3.18  yields 

E,  =  A±V,4>(x,y)exp(±jkz)  (3.29) 

3-4  Two-dimensional  Frequency  Limit 

In  solving  for  the  electric  field,  it  can  be  advantageous  to  remove  the  z  dependence  from 
the  electric  field.  Therefore,  let  us  consider  the  power  series  expansion  of  the  exp(±jkz)  terms  in 
Equation  3.29 

exp(±;*:r)  =  1  ±  (jkz)  +  ±  ■  +  •  ■  ■  (3.30) 

Thus,  if  kz  <C  1,  Equation  3.29  simplifies  to 

E,  wV„^(x,y)  (3.31) 

where  the  arbitrary  constant  (A±)  was  absorbed  into  <j>{x,y).  Equation  3.31  indicates  that,  for 
kz  I,  the  potential  in  the  tranverse  plane  is  essentially  constant  along  the  z  direction,  and  the 
subscript  t  notation  can  be  deleted  since  all  remaining  functions  are  transverse  functions.  Thus, 
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E  =  V<f>{x,y) 


(3.32) 


Since  the  wavenumber,  k,  is  related  to  frequency  by 


k  =  UJ^/JI€  = 


(3.33) 


there  is  a  specific  frequency  associated  with  the  condition  that  kz  =  which  limits  the  use  of 
the  approximation  of  Equation  3.31  (all  cell  layers  contain  nonmagnetic  material;  that  is  Hr  =  1). 


Since  c  = 


3  ■  10^°  cm/sec,  where  c  is  the  speed  of  light  in  a  vacviurn.  th*'  frequency  limit 


{uil)  imposed  due  to  the  length  of  the  fingers  (L)  is 


(3.34) 


For  the  interdigitated  electrode  in  the  CHEMFET,  L  ss  0.4  cm,  and  the  largest- valued  fr 
is  associated  with  layer  4  (mean  value  for  the  polyphthalocyanine  polymers  is  t,  «  650,  with  an 
overall  range  of  16-1300  [4:167]).  Thus,  for  this  case 


0.4v^ 


=  2.94  •  10®  rad/sec 


(3.35) 


which  translates  to  an  ordinary  frequency  limit,  //,,  associated  with  this  length  as 


fL  =  =  468  MHz 

27r 


(3.36) 


Thus,  for  frequencies  as  large  as  50  MHz  (/  <C  fi),  the  electric  field  in  the  CHEMFET’s  inter¬ 
digitated  electrode  structure  is  essentially  constant  along  the  z  direction,  thereby  transforming  the 
three-dimensional  problem  into  a  two-dimensional  problem  in  x  and  y. 


3.5  Laplace’s  Equation  in  the  Transverse  Plane 


In  the  transverse  plane,  the  electric  field  must  satisfy 


V,  •  E,  = 

f 
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where  p  is  the  charge  density.  In  a  charge  free  region  away  from  the  conductors,  this  equation 
simplifies  to 

V,  E,  =0  (3.38) 

The  subscript  t  notation  can  be  dropped,  since  the  only  variables  remaining  are  in  the  transverse, 
x-y,  plane  and  Equation  3.32  can  be  substituted  into  3.38  to  obtain  Laplace’s  equation. 

V^4,{x,y)  =  0  (3.39) 

3.6  Boundary  Conditions  in  the  Transverse  Plane 

In  the  transverse  plane,  the  voltage  must  be  continuous  at  the  boundary  between  two  layers 
(for  example  between  layers  1  and  2),  and  the  normal  component  of  the  E  field  must  satisfy 

fi-  El  -  ^e2  -  Ej  =0  (3.40) 

where  n  is  a  unit  vector  normal  to  the  boundary,  €„  is  the  permittivity  of  layer  n,  j  =  \/^,  (r„  is 
the  conductivity  of  layer  n,  uj  is  the  angular  (radian)  excitation  frequency,  and  E^  is  the  electric 
field  in  layer  n.  For  the  transverse  plane,  y  is  normal  to  the  layer  boundaries,  so  Equation  3.40 
reduces  to 

(<.-P)£„  =  (<,-;2)£„  ,3.4,) 

where  Eny  is  the  y  component  of  the  electric  field  in  layer  n. 

Since  this  relationship  is  frequency  dependent,  the  real  part  of  the  equation  will  dominate 
at  large  frequencies,  while  the  imaginary  part  will  dominate  at  low  frequencies.  The  high  fre¬ 
quency  behavior  of  Equation  3.41  can  be  demonstrated  by  considering  the  limit  as  w  becomes  large 
(lj  — ►  oo);  then,  Equation  3.41  reduces  to 

CiEiy  =  (2p7y  (3-42) 
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This  result  is  identical  to  the  case  with  no  conductivity  present  on  either  side  of  the  boundary.  By 
applying  this  high  frequency  condition  at  each  boundary,  the  cell  will  appear  to  contain  only  pure 
dielectric  material  in  each  layer  and  we  can  solve  for  the  potential  distribution  and,  eventually,  the 
capacitance  for  this  ideal  case.  Thus,  the  high  frequency  boundary  conditions  should  lead  to  the 
capacitive  elements  of  the  model. 

The  low  frequency  behavior  of  Equation  3.41  can  be  illustrated  by  multiplying  each  side  by 
u)  and  taking  the  limit  as  u;  — >  0  to  obtain,  for  example 


(TiEiy  —  a'2E2y 


(3.43) 


If  only  one  side  of  the  boundary  possesses  conductivity  (for  example,  layer  1),  the  above  equation 
simplifies  to  Eiy  =  0.  This  low  frequency  boundary  condition  permits  the  material  to  be  treated 
as  a  conductor.  Solving  for  the  potential  distribution  with  this  boundary  condition  will  eventually 
provide  the  resistance  of  the  material  in  this  ideal  case.  Since  the  sensor  coating  layer  (layer  4)  is 
surrounded  by  materials  with  <t  =  0  (air  and  insulator),  this  low  frequency  boundary  conditions 
can  be  enforced  on  layer  4  to  solve  for  the  potential  distribution  and,  eventually,  the  resistance 
through  this  layer  in  the  ideal  case. 

Within  the  cell,  the  conductivity  of  the  substrate  (layer  1)  is  much  greater  than  the  con¬ 
ductivity  of  the  sensor  coating  (layer  4).  This  situation  implies  that  the  low  frequency  boundary 
condition  can  apply  to  the  substrate  boundary  at  frequencies  where  the  high  frequency  boundary 
conditions  still  apply  to  the  sensor  coating  layer.  An  angular  frequency  (w,)  can  be  defined  to 
signify  the  point  where  the  substrate  transitions  from  a  conductor  to  a  dielectric  (the  point  where 
the  real  part  equals  the  imaginary  part  on  the  left  hand  side  of  Equation  3.41).  This  transition 
occurs  when 


w,  = 


(3.44) 
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By  restricting  the  model  to  frequencies  far  below  this  point,  the  substrate  can  be  treated  as  a  good 
conductor  to  further  simplify  the  model. 


For  example,  if  the  conductivity  of  the  substrate  is  0.2  15/cm  and  the  relative  permittivity  is 
11.9,  the  transition  frequency  is 


0.2 

“  11.9e„ 


(3.45) 


where  =  10~‘ Y36ff  F/cm  is  the  permittivity  of  free  space.  Thus,  the  radian  transition  frequency 


w,  m  1.9-  10^^  rad/sec 


(3.46) 


The  corresponding  transition  frequency  (/,)  is 

/,  =  =  30  GHz  (3.47) 

This  parameter  represents  the  frequency  where  the  substrate  layer  no  longer  acts  as  a  good  con¬ 
ductor.  Therefore,  if  the  operating  frequency  is  much  less  than  /, ,  the  substrate  remains  a  good 
conductor  (note:  Kraus  [13:446]  defines  a  good  conductor  as  lOOwe  <  cr).  For  example,  an  inter- 
digitated  electrode  structure  0.4  cm  long,  having  a  substrate  with  cri  =  0.2  U/cm  and  Cr,  =  11.9, 
will  have  /,  s;  100  /l.  Thus,  the  substrate  behaves  as  a  good  conductor  for  the  frequency  range 
(0  <  /  <  //,).  By  restricting  the  model  of  the  interdigitated  electrode  structure  to  frequencies 
below  the  transition  frequency  of  the  substrate,  the  ground  plane  appears  to  be  just  below  layer  2 
and  the  model  for  the  interdigital  sensor  structure  simplifies  to  that  shown  in  Figure  3.2. 

If  a  semiconductor  substrate  (silicon)  is  used,  a  depletion  layer  capacitance  must  also  be 
considered  [16:192-194].  This  capacitance  is  dependent  on  the  voltage  bias  to  the  metal-insulacor- 
semiconductor  contact  and  adds  a  significant  complexity  to  the  model  calculations.  This  additional 
capacitance  introduces  a  variable  model  element  (variable  depending  on  the  potential  of  each  elec¬ 
trode),  which  is  not  considered  in  this  thesis.  The  simplified  model  depicted  in  Figure  3.2  can  still 
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be  used  for  a  semiconductor  substrate,  if  a  suitable  ground  plane  is  introduced  between  layers  1 
and  2. 


Figure  3.2.  Simplified  model  for  the  interdigitated  electrode  sensor  structure 


The  sensor  coating  (layer  4)  has  a  much  smaller  conductivity  value  than  the  substrate.  There¬ 
fore,  throughout  the  operational  frequency  range  considered  herein  (0  <  f  <  fi  C  /»),  both  small 
and  large  w  behaviors  will  be  observed  at  the  boundaries  of  layer  4.  As  discussed  for  the  substrate, 
there  is  a  characteristic  frequency  associated  with  the  transition  of  layer  4  from  a  conductor  to  a 
lossy  dielectric.  Using  the  condition  where  <j^  =  WpfrA^o,  a  frequency  (wp)  can  be  defined  as  the 
angular  frequency  where  the  polymer  transitions  from  a  conductor  to  a  lossy  dielectric. 

For  the  sensor  coating  used  in  the  CHEMFET,  (Ta  =  10~*  U/cm  and  (a  =  650eo,  and 

0-4  10“^ 

u)p  =  —  =  -  =  174  rad/sec  (3.48) 

f4  650fo  '  ^  ' 

The  corresponding  ordinary  frequency  in  hertz  (/p)  is 
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(3.49) 


/p  =  ^  =  27.7  Hz 

Thus,  for  higher  frequencies  (fp  f  ft  ft),  the  sensor  coating  is  a  lossy  dielectric.  We 
can  solve  for  the  potential  and  the  capacitive  elements  of  Figure  3.2  for  the  ideal  case  of  a  pure 
dielectric  material  in  each  layer  by  enforcing  the  following  high  frequency  boundary  conditions  on 
the  potential,  according  to  the  development  of  Equations  3.41  and  3.42; 


—  ^zEzy 

(3.50) 

—  i^EAy 

(3.51) 

i^E^y 

—  ^sEzy 

(3.52) 

with  the  ground  plane  appearing  below  layer  2. 

Similarly,  for  low  frequencies  (0  <  /  fp),  we  can  solve  for  the  potential  in  the  conductive 
layer  4  and  model  it  by  a  lumped  resistive  element,  by  enforcing  the  following  low  frequency 
boundary  conditions  on  the  potential  as  in  Equation  3.43: 

£4,  =0  (3.53) 

3.7  Analysts  Summary 

In  summary,  to  find  the  electrical  parameters  of  interest,  we  must  solve  for  the  potential 
distribution  in  the  transverse  plane.  This  potential  function  is  periodic  and  incorporates  just  four 
material  layers  (the  ground  plane  appears  above  layer  1  since  the  substrate  is  a  good  conductor  for 
the  frequency  range  of  interest).  The  resulting  problem  diagram  is  shown  in  Figure  3.3. 

In  this  diagram,  d  is  the  insulator  thickness,  h  is  the  electrode  height,  /  is  the  combined  height 
of  the  electrode  and  sensor  coating  layer,  while  s  is  the  electrode  separation.  These  quantities 
are  normalized  to  the  width  of  the  electrode  and  are,  hence,  nondimensional  quantities.  Vi  and 
V2  represent  the  potential  on  the  left  and  right  electrode,  respectively.  Within  each  layer,  the 
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potential  must  satisfy  Laplace’s  equation.  In  addition,  the  potential  above  layer  5  must  be  bounded 
at  infinity.  Using  the  boundary  data  shown  in  the  diagram,  the  potential  must  then  satisfy  the 
following  conditions: 


= 

0 

for  k  =  2, 3, 4, 5 

(3.54) 

d4)k,s+\ 

ax  2  ’ 

= 

0 

for  ib  =  2, 4,  5 

(3.55) 

4>7{x,  -d) 

= 

0 

for  |x|  < 

(3.56) 

lim  |<i5(x,y)| 

y  — *00 

< 

oo 

for  |x|  < 

(3.57) 

o!>2(^,0) 

= 

U, 

for  <  X  <  -f 

(3.58) 

h) 

= 

U, 

for  <  X  <  -| 

(3.59) 

(/>2(x,  0) 

U2 

for  1  <  X  <  ^ 

(3.60) 
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h)  = 

Vn  for  i  <  X  < 

(3.61) 

03(-^. 

■  y)  = 

V'\  for  0  <  y  <  /i 

(3.62) 

03(^. 

y)  = 

V2  for  0  <  y  <  /i 

(3.63) 

In  addition,  the  voltage  must  satisfy  the  boundary  conditions  at  the  interfaces  between  the  layers. 
First,  the  potential  is  continuous  across  layers: 


0) 

=  <i>3(x,0) 

for 

ki<  § 

(3.64) 

<i'3(-c,  h) 

=  <t>4{x,h) 

for 

ki<  1 

(3.65) 

-  4>5{x,1) 

for 

|x|< 

(3.66) 

Also,  the  electric  field  normal  to  each  interface  must  satisfy  equations  analogous  to  Equation  3.41. 

To  determine  the  capacitive  elements,  we  solve  for  the  potential  distribution  using  the  high 
frequency  boundary  conditions  of  the  normal  component  of  the  electric  field: 


(7E‘2y 

for 

|x|  <  1  and  y  =  0 

(3.67) 

(sEsy 

for 

II 

a 

V 

(3.68) 

^^E^y 

II 

for 

(i)  <  and  y  =  1 

(3.69) 

Since  Ey  =  —d<t)ldy  in  each  of  the  layers,  these  equations  yield  the  following  conditions  on  the 
potential  in  each  of  the  layers: 


^(x.O)  = 

94)3, 

-tr-C'T.O 

dy 

for  |x|  <  ^ 

(3.70) 

84)2 ,  , , 

d4>4  , , 

for  |x|  <  § 

dy 

dy 

(3.71) 

84>4 1  ,, 

= 

d4>5,  ,, 

for  |x|  <  ^ 

(3.72) 
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To  determine  the  resistive  element,  we  solve  only  for  the  potential  distribution  in  layer  4,  using 
the  low  frequency  boundary  conditions  of  the  normal  component  of  the  electric  field,  Equation  3.43. 
From  Figure  3.3,  with  Vi  =  +1/2  and  V2  —  —lf2,  the  potential  distribution  must  satisfy  the 
following  conditions: 


=  0 

(3.73) 

4>4(x,h)  = 

for  <  X  <  -| 

(3.74) 

<t>A{x,h)  =  i 

for  f  <  X  <  iiT 

(3.75) 

904  f  s+l  \ 

2  ’V  =  " 

for  /i  <  y  <  / 

(3.76) 

The  low  frequency  boundary  conditions  of  the  normal  component  of  the  electric  field  are 

Ey  =  0  for  |zl  <  and 

II 

(3.77) 

Ey  =  0  for  jij  <  and  y  =  1 

(3.78) 

Since  Ey  =  —d4>\ldy,  these  equations  transform  into  the  following  conditions  on  the  potential  in 
layer  4: 


904 

^(x,h)  =  0  for|x|<i 

(3.79) 

^(x,/)  =  0  for  |x|  < 

(3.80) 
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IV.  Capacitance  Solution 


To  determine  the  capacitive  elements  of  the  model,  we  solve  Laplace’s  equation  for  the  po¬ 
tential  distribution  within  the  periodic  cell,  using  the  boundary  conditions  corresponding  to  the 
large  frequency  limit.  First,  the  potential  functions  of  each  layer  are  found  which  satisfy  the  partial 
differential  equation  and  boundary  conditions  at  x  =  ±(s  +  l)/2  or  r  =  ±s/2,  as  appropriate. 
These  functions  contain  constants,  which  are  determined  by  enforcing  the  boundary  conditions  at 
the  interfaces  between  the  layers.  The  electric  field  is  then  determined  from  the  gradient  of  the 
potential,  and  the  charge  is  found  using  Gauss’  law.  The  capacitance  values  are  found  directly  from 
the  charge  using  Q  =  CV ,  where  Q  is  charge,  C  is  the  capacitance  of  interest,  and  V  is  the  known 
potential.  A  diagram  depicting  the  simplified  periodic  cell  and  the  notation  is  shown  in  Figure  4.1. 


Figure  4.1.  Problem  diagram  for  capacitive  solution 
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4-1  Potential  in  Each  Layer 


The  potential  function  in  each  region  was  determined  using  the  separation  of  variables  tech¬ 
nique  with  Laplace’s  equation  in  the  transverse  (x-y)  plane.  Eigenfunctions  of  the  independent 
variable  x  are  needed  to  match  the  potential  and  normal  electric  field  across  layer  interfaces.  Pre¬ 
sented  here  is  a  sumary  of  the  solution,  which  is  detailed  in  Appendix  A. 

4.1.1  Potential  Functions  Using  the  separation  of  variables  technique,  we  obtain  the  po¬ 
tential  4>k  in  layer  k  Lot  k  =  2,  3, 4,  5  as  functions  of  the  impressed  voltages,  Vi  and  V2,  and  all  cell 
parameters.  These  potential  functions  are; 


<i>:>{x,y) 

4>A{x,y) 

4>5(x,y) 


00 

ao(y  -t-  d)  -I-  ^  a„  sinh[A2(j/  +  d)]i!>n{x) 

n  =  l 

00  je 

bn  cosh 

n  =  l 

00 

n  =  1 
00 

do  +  d„  exp(-A5y)V'®  (f) 

n  =  l 


-1-  Cn  sinh 


"'2 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


where 


At 

A3 

v^i(x) 

Kiy) 


HTT 

S  +  1 

riTT 

s 


cos 


flTT 


for  A:  =  2,  4,  5 


X  1 

7+1  2 


sin 

'U2-  Ui 


1  *  1 

"'17  +  5 


X  + 


for  it  =  2, 4, 5 


V  2  +  Pi 


exp(-A4/)  |cosh[A4(/  -  j/)]  -t  ^  sinh[A4(/  -  j/)]| 


(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
(4.10) 


and  On,  bn,  Cn,  and  dn  are  constants  used  to  satisfy  the  boundary  conditions.  The  superscript  in 
tniy)  denotes  the  layer  of  interest,  instead  of  the  normal  “power”  interpretation. 
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4- 1.2  Applying  Boundary  Conditions  to  Determine  the  Coefficients.  To  determine  the  un¬ 
known  constants,  we  have  the  boundary  conditions  to  enforce  at  the  interface  between  layers  2  and 
3  and  that  between  layers  3  and  4.  The  potential  is  continuous  across  each  layer  (continuity  con¬ 


dition),  while  the  y  component  of  the  electric  field  must  satisfy  the  capacitive  boundary  conditions 
imposed  by  Equations  3,70,  3.71,  and  3.72  (jump  condition).  Therefore,  we  use  the  orthogonal 
properties  oi  ihe  irigonomeiric  functions  to  “sift"  certain  coefficients  out  of  the  infinite  series.  This 
will  allow  us  to  develop  a  system  of  equations  which  we  can  solve  to  determine  the  four  scis  of 
coefficients. 

For  the  continuity  condition,  we  sift  with  the  outer  (layers  2  and  4)  eigenfunction  since  the 
boundary  data  is  known  across  the  entire  width  of  the  periodic  cell.  However,  we  sift  with  the  inner 
(layer  3)  eigenfunction  for  the  jump  condition,  since  it  is  restricted  to  |2;|  <  From  the  solution 
detailed  in  Appendix  A,  the  following  equations  must  be  solved  to  find  the  coefficients  a„,  6„,  c„, 
and  dn'. 


aod{s  -b  1)  —  “  (“0"]  +  <^o(s  +1)  —  (^  +  1)  (T)  +  V2) 


nsl 


aod(s-fl)  +  ^ — Cnip' [1  —  (  — 1)"]  —  <fo(s  +  1)  =  0 


n=:l 


-  00 


n  =  l 


t 

+  2  C.niL  - 


=  0 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


n  =  l 


lif. 

mTT 

111 

TTlTr 


[1  -  (-1)"']  ao  +  Y.  ^  =  0  (4.15) 

,  s  +  1  “  s  +  1 

n=l  n=l 

00  00 

[1  -  (-l)'^]ao  +  Y  -  Y  =  0  (4,16) 


where  m  is  a  positive  integer  and 


©L  =  sinh 


(mird 

tttJ 


(4.17) 
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r'  =  cosh 


(mird 


=  exp  (-Am/) 

KW 

=  exp(-A„/) 
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^  n 

.  ,  f nvh^ 
—  smh  - 

V  2s  , 

A„ 

nn 
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(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 
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(4.24) 
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for  n  = 

»+i 


-t=i”(jrri) 


(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 


4.2  Charge  Calculations 


To  find  the  charge  on  either  microstrip,  we  integrate  the  normal  electric  flux  density  around 
a  closed  contour  enclosing  the  microstrip.  The  result  of  the  charge  calculations  is  simply  presented 
here,  while  the  calculations  are  detailed  in  Appendix  B.  The  charge  on  either  microstrip  is 
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^C2ao(s  +  1)  +  €2  X]  (^)  ^  “  ^^2) 

n  =  l 

00  /  L\  ^ 

-2(3  cos  ^  y)  (^)  +  H  (  y) 

n=l  ^  '  n=l 


(4.30) 
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4-3  Fxaluation  of  the  Capacitive  Elements 


To  solve  for  the  capacitive  elements,  we  utilize  Weeks’  method  [17]  to  solve  for  the  elements 
one  at  a  time.  The  capacitance  to  ground  element,  Cj,  is  found  by  setting  both  electrodes  to 
the  same  potential  (Vi  =  V2).  The  charge  on  either  conductor  can  then  only  be  associated  with  a 
capacitance  to  ground,  since  an  electric  flux  cannot  exist  between  the  two  conductors  when  they  are 
at  the  same  potential.  Therefore,  if  we  solve  for  the  potential  and  charge  for  the  case  of  Vi  =  V2  =  1 
volt,  then 


Cg  -  <3(Vi=Vj=l) 


(4.31) 


where  Q(Vi=V'j=i)  is  the  charge  for  the  case  of  Vi  —  =  I  volt. 

To  solve  for  Cg,  one  of  the  electrodes  is  grounded  (for  example,  Tj  =  1  volt  and  V2  =  0 
volts).  Now  the  total  capacitance  to  ground  (Ct),  from  the  ungrounded  electrode,  is  the  parallel 
combination  of  Cg  and  Cg-  Since  capacitances  in  parallel  add,  then 

Ce  =  Ct-C, 

=  <3(V',  =  l,Vj=0)  -  Q(J^,  =  V,  =  1) 

where  Q(Vi  =  i,V3=o)  's  the  charge  for  the  case  of  Vj  =  1  volt  and  V2  =  0  volts. 

This  method  requires  the  potential  to  be  solved  for  two  separate  cases.  Therefore,  it  can  be 
advantageous  to  use  symmetry  within  the  periodic  cf'.l  to  reduce  the  computational  requirements 
and  solve  for  the  potential  only  once.  Hence,  it  is  beneficial  to  consider  symmetry  within  the 
periodic  cell. 

4-3.1  Cell  Symmetry  and  the  Eigenfunctions.  To  investigate  possible  symmetry  within  the 
periodic  cell,  it  is  necessary  to  further  analyze  the  eigenfunctions  in  each  layer  of  the  cell.  From 
Equation  A.  16,  the  eigenfunctions  in  layers  2,  4,  and  5  are 


(4.32) 

(4.33) 
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'l'n(x)  =  yl„  cos  HT  (  — ^  +  i 

\  s  +  1  2 


(4.34) 


Since 

cos(a  +  /3)  =  cos(a)  cos{l3)  —  sin(Q')  sin(/?) 
we  may  write 

'I'r.W 

where 


=  An  cos 


nnx 


=  'lr‘ 


s  +  1 
(x)  - 


A  ■  (  ””  ^ 

I  —  1  —  An  Sin  I  -  )  sin  1  — 

\  2  }  \s+lj  W/ 


(4.35) 


(4.36) 

(4.37) 


(4.38) 

(4.39) 


Hence,  the  original  eigenfunction  has  components  which  have  odd  and  even  symmetry  with  respect 
to  X  =  0.  If  the  applied  voltages  are  equal,  then  the  potential  must  be  symmetric  with  respect  to 
X  =  0,  so  that 


'I'„(x)  =  '^''''"(x)  (4,40) 

and 

(4.41) 

Thus,  all  the  coefficients  of  the  odd-indexed  terms  must  be  zero  when  the  potential  has  even 
symmetry  with  respect  to  x  =  0,  due  to  the  factor  sin(n)r/2)  in  Similarly  with  odd  symmetry 

=  -V2) 
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and 


=:  0  (4.43) 

Hence,  all  the  coefficients  of  the  even-indexed  terms  must  be  zero  when  the  potential  has  odd 
symmetry  with  respect  to  x  =  0,  due  to  the  factor  cos(mr/2)  in 

For  layer  3,  from  Equation  A. 50,  the  eigenfunction  is 

$^(x)  =  fl„sin  "’f  5  j  (4.44) 

Now,  employing  the  trigonometric  identify 


sin(a  -(-/?)=:  sin(at)  cos(/?)  +  cos(a)  sin(/?)  (4.45) 

we  have 

rT.  ^  ^  •  /n7rx\  ^nir\  „  /nKX\  .  fnv\ 

^^(x)  =  S„sin  cos  j -f  fl„cos(^-yj  sin  j  (4.46) 

_  ^  ^^."'"(x)  (4.47) 

where 

(~)  (^) 

'l’^'"'"(x)  =  B„cos(^)sin(^)  (4.49) 

With  even  symmetry  (V^i  =  V2) 

'l'^(x)  =  '^^'"(x)  (4.50) 


and 

=  0  (4.51) 
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Thus,  all  the  coefficients  of  the  even-indexed  terms  in  layer  3  must  be  zero  for  even  symmetry.  For 
odd  symmetry  (Vi  =  —  V2) 

^3(^)  ^  (4^52) 


and 


=  0  (4.53) 

Hence,  all  the  coefficients  of  the  odd-endexed  terms  in  layer  3  must  be  zero  for  odd  symmetry. 

4-3.2  Symmetry  and  the  Charge  Function.  We  can  use  tbe  result'  of  the  previous  section 
to  simplify  the  charge  expression.  From  Equation  4.30 

Q  =  ^f2ao(s+ +  (7^) 

n  =  l  ' 

-2£3  cos  (y )  6„  sinh  -F  £4  £  (y) 

n=l  '  '  n=l 

=  Q1 -FQ2  4-Q3-|-g4-|-Q5  (4.55) 

where  each  term  in  order  was  assigned  a  corresponding  label.  With  even  symmetry  (Kj  =  V2),  the 
odd  terms  in  layers  2  and  4  are  zero  (a„  =  d„  =  0  for  odd  n)  and  the  even  terms  in  layer  3  are  zero 
(hr,  =  c„  =  0  for  even  n).  Hence, 

Q2  =  Q3  =  Q4  =  Q5  =  0  (4.56) 

and 

Q— =  Q1  =  ir2ao(s  +  l)  (4.57) 

For  odd  symmetry  (V)  =  —V2),  the  even  terms  in  layer  2  are  zero;  thus,  uq  =  0  and 

=  Q2 -F  Q3  +  Q4 -f  Q5  (4.58) 
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These  symmetry  considerations  are  important  because  they  simplify  the  calculations  required 
for  the  capacitive  elements  of  the  model.  To  determine  the  Cg  element,  both  electrodes  are  fixed 
to  the  same  potential  (for  example  Vj  =  V2  =  1  volt),  while  for  the  Ce  element,  one  electrode  is 
grounded  (for  example  V)  =  1  volt  with  V2  =  0).  However,  the  non-symmetric  case  (l/i  =  1  and 
Vo  =  0)  is  equivalent  to  the  superposition  of  even  symmetry  (Tj  =  Vo  =  .5)  and  odd  symmetry 
(Vi  =  — V2  =  .5).  Thus,  using  cell  symmetry,  we  have  a  method  to  extract  the  even  symmetry 
charge  using  the  solution  to  the  potential  in  this  non-symetric  case. 

Therefore,  we  only  have  to  solve  for  the  potential  coefficients  once  (in  the  non-symmetric 
case)  to  determine  both  capacitive  elements.  Without  these  symmetry  arguments,  the  potential 
coefficients  would  have  to  be  determined  separately  for  each  case. 

4-3.3  Capacitance  to  Ground:  Cg.  With  even  symmetry  in  the  transverse  plane,  both  elec¬ 
trodes  have  the  same  potential,  and  the  charge  on  either  conductor  will  only  be  associated  with  a 
capacitance  to  ground.  Hence  this  capacitance  to  ground  Cg  can  be  found  from  the  relationship 

Cg  =  (4.59) 

where  Q'"'"  is  the  charge  from  the  even  symmetry  case  and  V  is  the  potential  on  the  electrode.  The 
potential  coefficients  are  calculated  using  V(  =  1  and  V^  =  0.  Since  this  situation  is  equivalent  to  the 
superposition  of  an  even  symmetry  case  (V'l  =  V2  =  .5)  and  an  odd  symmetry  case  (Vi  =  —V^  =  .5), 
we  use  T  =  .5  in  the  formula  corresponding  to  the  value  in  the  even  symmetry  crise.  Thus 

Cj  =  2Q''''"  =  e2ao(s+ll  (4.60) 

4-3.4  Capacitance  Between  Electrodes:  Ce-  With  Vi  =  I  and  V2  =  0,  the  total  capacitance 
to  ground  will  be  the  parallel  combination  of  Ce  and  Cg.  Since  capacitances  in  parallel  add,  the 
total  capacitance,  Ct,  is 

CT  =  CE+Cg  =  ^  (4.61) 
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where  Q  is  the  total  charge  and  V  is  the  potential  on  the  electrode.  Hence,  with  V'  =  1 


Ce  =  Ct  -  Cg=Q-  =  Q'"'”'  -  Q*"*" 


or 


(4.62) 


Ce 


€2  ^  a„  sin  r'  +  ^  (Vi  -  Va)  -  2f3  £  cos  (^) 

n  =  l  n  =  l  '  ■' 

oo  - 

+  f4  (^) 

n=2l 


4-4  Estimating  the  Potential  Coefficients 


By  truncating  the  infinite  series  of  Equations  4.11-4.16  at  N  terms  (where  N  is  an  integer), 
we  can  estimate  the  first  N  terms  of  each  set  of  potential  coefficients.  Therefore  we  can  deter¬ 
mine  ao,  Ui, . . . ,  as-i,  61,62,  •  •  -  iV,  ci,C2, ..  .,c/v,  and  do.di, . . . ,  ds-i  by  solving  a  system  of  iN 
Equations.  The  equivalent  matrix  equation  formulation  is 


" 

• 

- 

• 

61mn 

0 

dim 

dn 

Fm 

^  Im 

0 

clmn 

-dim 

6„ 

0 

b2m 

0 

d2mn 

0 

0 

c2m 

-d2mn 

d„ 

0 

where 


al 


m 


61 


mn 


clmn 


d{s  -1-1)  for  m  =  0 
form  5^0 

V. 

[  -^’^n[l-(-ir]  form  =  0 
I  -2«I>=/^„  form  5^0 

I  form  =  0 

2't'/3„  form  5^0 


(4.64) 


(4.65) 


(4.66) 


(4.67) 
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I 

I 


fi\m 

^  1 

5-1-1  for 

m  =  0 

1 

Kih)^  for 

t 

3 

-H. 

o 

(4.68) 

^2j7in 

' 

^[i-(-in 

for  n  =  0 

=  < 

l2lLlrc  t3 
s  +  i  ^  n  ^nm 

for  n  ^  0 

(4.69) 

b2m 

(4.70) 

c2fji 

=  - 

-famTr*^^ 

(4.71) 

_  f^nTT  4 

'  «  +  i’'" 

(4.72) 

(s  -|-  l)(Vi  -f-  Vo) 

for  m  =  0 

=  1 

2Tm 

for  m  ^  0 

(4.73) 
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V.  Resistive  Solution 


T'l  determine  the  resistive  element  of  the  model,  we  solve  Laplace’s  equation  for  the  potential 
distribution  within  the  conductive  layer  4,  using  the  boundary  conditions  associated  with  the  low 
frequency  limit.  First,  the  layer  is  divided  into  three  regions  to  match  the  boundary  data.  The 
potential  functions  in  each  region  are  then  found,  which  satisfy  the  partial  differential  equation  and 
the  boundary  conditions  at  y  =  /»  and  y  =  1.  These  functions  contain  unknown  constants  that  are 
determined  by  enforcing  boundary  conditions  at  the  interfaces  between  the  regions.  The  electric 
field  (E)  is  determined  from  the  gradient  of  the  potential  and  the  current  density  (J)  is  found 
using  J  =  crE,  where  cr  is  the  conductivity  of  the  layer.  The  current  through  the  layer  is  found  by 
integrating  the  current  density  over  a  suitable  contour,  and  the  resistance  is  then  determined  using 
Ohm’s  law.  The  problem  diagram  is  shown  in  Figure  5.1,  in  which  the  x  axis  has  been  translated 
by  -h  to  simplify  the  analysis.  Thus,  0  <  y  <  p,  where  p  =  I  —  h. 


A  y 

p  1 

"  ^  1 

]  region  a 

1 

1 

region  6 

region  c  i 

1 

1 

!  =  .5 

l' 

II 

0 - - 

-(s+l)/2  -s/2  0  s/2  (s+l)/2 


Figure  5.1.  Diagram  of  resistive  layer 


•5. 1  Resistive  Potential  Functions 

The  resistive  layer  is  divided  into  3  regions  (see  Figure  5.1).  Region  a  is  the  left  region 
which  is  bounded  by  — (s  l)/2  <  x  <  —sf2.  Region  6  is  the  center  region  which  is  bounded  by 

i 


—s/2  <  X  <  s/2.  Region  c  is  the  right  region  which  is  bounded  by  s/2  <  x  <  (s  +  l)/2.  In  the 
resistive  solution,  we  concentrate  on  the  case  of  odd  symmetry  with  respect  to  x  =  0  (Vi  =  —  Va  = 
1/2),  since  the  even  symmetry  case  (V'l  =  Vj)  does  not  produce  any  electric  flux  between  the  two 
electrodes.  Using  the  resistive  boundary  conditions  and  the  separation  of  variables  approach,  a 
series  form  of  the  voltage  potential  is  obtained,  with  eigenfunctions  in  the  independent  variable 
y,  consisting  of  trigonometric  functions,  while  the  corresponding  functions  in  x  are  hyperbolic 
functions.  For  simplification,  the  results  are  summarized  here,  while  the  complete  solution  is 
detailed  in  Appendix  C  (<!)*  denotes  the  potential  in  region  k,  with  k  =  a,  b,  or  c).  The  potentials 
in  the  res'stive  solution  are 


4>aix,y) 

<t>bix,y) 


1 

2  +  e„  cosh 

nsl 

00 

fo  +  '^gn  sinh 

n  =  l 


sin 


(2n  -  \)iry 

2p 


‘i>c{x,y)  =  -</>a(-x,y) 


(5.1) 


(5.2) 

(5.3) 


where  e„,  /o ,  and  g„  are  constants  used  to  satisfy  the  continuity  of  the  potential  and  its  x  derivative, 
across  the  boundaries  of  the  regions. 


5.2  Applying  Boundary  Data  Between  Regions 

The  constants  can  be  determined  by  enforcing  the  boundary  conditions  between  the  regions. 
Because  of  symmetry  with  respect  to  x  =  0,  either  the  left  boundary  (regions  a  and  6)  or  the  right 
boundary  (regions  b  and  c)  produce  the  same  system  of  equations  to  solve.  At  the  boundary,  the 
potential  is  continuous  across  the  boundary  (continuity  condition)  and  the  x  component  of  the 
electric  field  is  also  continuous  (jump  condition)  since  the  conductivity  is  identical  on  both  sides  of 
the  boundary.  This  situation  leads  to  two  boundary  equations  with  infinite  series  on  each  side  of  the 
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equations.  However,  orthogonality  of  the  cosine  functions  can  be  used  to  reduce  the  two  equations 
to  only  an  infinite  series  involving  the  e„  coefficients.  The  resulting  equations  from  Appendix  C 
are 


n  =  1 

4(2n-l)  2(2n-l)^ 

x[(2n- l)2-dm2]  "  7r((2n-l)2-4m2]  ” 


pr  ^ _ 1- 
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-1 
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where  m  is  a  positive  integer  and 
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[(2n- 

4p 
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=  sinh 
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f  rriTTs'' 

1  2p  > 

'Hm 

=  sinh 

( rmcs\ 
\  2p 

(.5.4) 

(5.5) 


(5.6) 


(5.7) 

(5.8) 

(5.9) 


The  superscript  r  notation  denotes  quantities  used  in  the  resistive  solution.  Once  the  e„’s  are 
determined,  the  gm's  and  fa  can  be  found  from  the  relationships 


fo 


9m 


1 
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1  r  4(2n  -  1) 

"T[(2n-1)3-W] 


(5.10) 
(5  11) 


5.3  Restsiance  Calculation 

To  determine  the  resistance  of  layer  4,  we  calculate  the  current  flowing  through  layer  4.  From 
Ohm’s  law 

V  -  IR  (5.12) 

where  V  is  the  voltage  potential.  I  is  the  current,  and  R  is  the  resistance.  Equivalently, 
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G=[/V  (5.13) 

where  G  =  \/R,  and  G  is  the  conductance  of  layer  4.  The  current  density  (J)  is  related  to  the 
electric  field  (E)  by 

J  =  (rE  (5.14) 


where  a  is  the  conductivity  within  layer  4.  The  total  current  can  be  determined  by  integrating  the 
current  density  over  a  suitable  contour,  which  encompasses  all  the  current  density.  Referring  to 
Figure  5.1,  the  current  density  in  the  x  direction  (Jj)  is  integrated  along  the  y  axis  from  y  =  0  to 
y  =  p.  From  Equation  5.2,  the  potential  in  region  6  is 


OO 

y)  =  /o-r  +  X]  sinh 

n  =  l 


Using  the  relationship 
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dz 


(5.16) 


we  may  write 


y  ^  TITT 
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^  T> 


nTl  P 


(5.17) 


,‘\t  X  =  0,  this  reduces  to 


Ei(G,y)  -  -h (5.18) 

P  VP/ 

Using  Equation  5.18  in  Equation  5.14  yields  the  current  density  along  the  y  axis,  where  x  =  0. 
Thus, 
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The  current  through  layer  4  is,  formally, 


(5.20) 

Using  Equation  5.20  in  5.13,  with  V  =  1  volt,  yields 

G=-<TfoP  (5.21) 

Equivalently, 

R=(-afopr'  (5.22) 

5.4  Estimating  the  Coefficients 

With  the  equations  presented  in  the  section  5.2,  the  unknown  constants  can  now  be  estimated 
by  truncating  the  infinite  series  in  Equations  5.4  and  5.5  at  an  integer  N .  Converting  the  equations 
to  matrix  form  yields 

[.V/mn][en]  =  [-Fm]  (5.23) 
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for  rn  ^  0 


(5.24) 


—  1  for  m  =  0 
0  for  m  0 


(5,25) 


Hence  we  can  approximate  ei ,  ct,  ■ . . ,  cat  by  solving  N  equations  (m  =  0,  1 , 2, . . . ,  yV  —  1),  and 
thus  obtain  an  NxN  matrix.  From  the  e„  terms,  we  can  also  estimate  /o  and  the  gm  terms  using 
Equations  5.10  and  5.11. 
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VI.  Numerical  Results 


Up  to  this  point,  the  analysis  and  solutions  have  been  general  in  nature,  with  dimensional 
and  material  values  of  the  periodic  cell  unspecified.  In  this  chapter,  we  use  the  expressions  de¬ 
veloped  in  previous  chapters  with  particular  values  for  the  structure’s  dimensions  and  material 
properties.  Specifically,  we  focus  on  cell  dimensions  and  material  aspects  which  are  representative 
of  the  chemically-sensitive  field-effect  transistor  (CHEMFET)  [18].  Within  this  chapter,  section 

6.1  concentrates  on  the  results  of  the  capacitive  solution,  while  section  6.2  focuses  on  the  resistive 
solution. 

In  sections  6.1.1  and  6.1.2,  the  capacitive  solution  is  illustrated  by  investigating  the  potential 
along  the  interfaces  between  layers  and  by  observing  constant  voltage  contour  plots  using  special 
test  data  for  the  cell  parameters.  In  section  6.1.3,  graphs  of  the  capacitive  elements  are  presented 
using  a  baseline  set  of  cell  parameters,  while  varying  individual  parameters  over  a  specific  range. 

In  section  6.2.1,  we  illustrate  the  resistive  solution  by  observing  the  potential  along  the  lower 
border  and  a  constant  voltage  contour  plot  created  using  special  test  values  for  the  layer  parameters. 
In  section  6.2.2,  graphs  of  the  resistive  elements  are  presented  using  the  baseline  values  for  the  cell, 
while  varying  individual  parameters  over  a  specific  range. 

6.1  Results  of  the  Capacitive  Solution 

A  FORTRAN  computer  program  wats  created  to  solve  the  truncated  linear  system  of  equations 
discussed  in  section  4.4,  and  thereby  approximate  the  coefficients  involved  in  the  capacitive  solution 
of  the  potential  function.  The  program  utilizes  the  IMSL  gaussian  elimination  FORTRAN  routine 
“LSARG”  [7:11]. 

6.1.1  Potential  Along  the  Boundaries.  By  truncating  the  series  at  A’  terms  (where  N  is 
an  integer),  we  are  able  to  estimate  the  first  N  coefficients.  For  the  capacitive  solution,  these 
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coefficients  are:  ao,  aj , . . . ,  ajv-i,  6i ,  62, . . . ,  6yv,  Ci,  C2, . . . ,  cjv,  and  do,  di , . . , ,  d^-i-  As  ;V  becomes 
larger,  the  overall  accuracy  of  the  final  sum  will  improve,  but  at  the  cost  of  increased  computer 
processing  time  to  solve  the  coefficients.  A  worthwhile  objective  is  to  use  the  fewest  number  of 
coefficients  (smallest  N)  and  yet  achieve  accurate  results. 

To  select  a  suitable  number  of  terms,  N ,  the  potentials  along  the  lower  boundary  (between 
layers  2  and  3)  and  the  upper  boundary  (between  layers  3  and  4)  are  investigated  with  N  =  10, 
20,  and  50.  Normalized  data  values  used  for  the  program  are  displayed  in  Table  6.1  (cr*  is  the 
relative  permittivity  or  dielectric  constant  in  layer  k,  and  =  frifo  where  to  is  the  permittivity 
of  free  space).  The  values  for  the  cell  parameters  were  selected  using  the  CHEMFET  data  as  a 
guide  and  normalizing  all  cell  dimension  with  respect  to  the  electrode  width.  The  results  at  the 
lower  boundary  are  displayed  in  Figure  6.1,  while  the  results  at  the  upper  boundary  are  shown  in 
Figure  6.2. 

Table  6.1.  Normalized  data  used  for  boundary  potentials  (with  no  coating) 


parameter 

d 

h 

1 

s 

fra 

^^4 

value 

0.1 

0.1 

Oil 

1 

3,9 

1 

1 

1 

With  N  =  50,  both  Figures  6.1  and  6.2  show  good  agreement  with  Equations  3.58,  3.59,  3.60, 
and  3.61  for  Vj  =  1  and  V2  =  0.  There  is  some  deviation  at  the  edge  of  the  left  electrode,  but  this 
is  expected  due  to  the  finite  number  of  terms  kept  in  the  Fourier  representation  of  the  potential 
function.  The  potentials  along  the  boundaries  are  non-linear  due  to  the  small  insulator  thickness  (d). 
Figure  6.3  displays  both  the  upper  and  lower  boundary  potential  with  various  insulator  thicknesses. 

With  Sensor  Coating.  Figure  6.4  shows  the  potential  function  along  the  lower 
and  upper  boundaries  with  various  insulator  thicknesses,  d,  using  the  data  in  Table  6.2.  Again, 
the  cell  parameter  values  were  selected  using  the  CHEMFET  structure  cis  a  guide  and  normalizing 
the  cell  dimensions  with  respect  to  the  electrode’s  width.  At  the  upper  boundary,  the  potential  is 
fairly  linear  from  one  electrode  to  the  other  (with  d  =  5).  The  potential  is  perturbed  from  this 
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(a)  Lower  Boundary  with  iV  =  10 


X  position 


(b)  Lower  Boundary  with  AT  =  20 


X  position 


(c)  Lower  Boundary  with  =  50 


Figure  6.1.  Potential  function  along  lower  boundary  with  N  =  10,  20,  and  50 
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(a)  Upper  Boundary  with  =  10 


X  position 

(b)  Upper  Boundary  with  N  =  20 


X  position 

(c)  Upper  Boundary  with  =  50 


X  position 


Figure  6.2.  Potential  function  along  upper  boundary  with  N  —  10,  20,  and  50 
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linear  behavior  as  the  insulator  is  made  thinner  (smaller  d).  At  the  lower  boundary  with  d  =  5,  the 
potential  is  still  nearly  linear,  but  the  linearity  is  much  more  perturbed  with  a  thinner  insulator 
layer  compared  to  the  upper  boundary. 

Comparing  Figures  6.3  and  6.4,  the  potential  is  nearly  linear  with  the  sensor  coating  above 
the  electrodes  as  compared  to  the  case  of  just  air  covering  the  electrodes. 


Table  6.2.  Normalized  data  used  for  boundary  potentials  (with  coating) 


parameter 

KB 

BK 

1 

S 

EBI 

mm 

mm 

E9II 

value 

0.1 

0.1 

0.11 

1 

3.9 

3.9 

500 

1 

6.1.2  Constant  Voltage  Contours.  A  FORTRAN  computer  program  was  created  to  calculate 
the  potential  in  the  cell,  using  the  voltage  coefficients  obtained  from  the  IMSL  routine  LSARG.  A 
100  by  100  point  grid  was  established  to  view  the  potential  in  the  cell,  with  a  spacing  of  (s  +  1)/100 
between  grid  points.  The  voltage  contour  plots  were  created  using  AFIT’s  METALIB  FORTRAN 
plotting  routines;  specifically,  the  contour  plotting  routine  “RCONTR”  [1]  was  utilized.  Constant 
voltage  contour  plots  were  created  for  five  cases.  The  model  parameters  of  the  five  cases  are  shown 
in  Table  6.3.  The  fundamental  cell  parameters  were  chosen  using  the  CHEMFET  structure  as 
a  guide;  however  specific  variations  were  selected  in  each  case  to  highlight  certain  aspects  of  the 
solution.  Here  again,  the  cell  dimensions  were  normalized  with  respect  to  the  electrode’s  width. 


Table  6.3.  Data  used  for  voltage  contour  plots 


case 

V^i 

Vz 

d 

mm 

1 

D 

mm 

mm 

A 

6.5 

1 

0 

1 

0.1 

0.11 

1 

1 

1 

1 

1 

B 

6.6 

1 

-1 

1 

0.1 

0.11 

1 

1 

1 

1 

1 

C 

6.7 

1 

-1 

1 

0.11 

1 

1 

100 

1 

1 

D 

6.8 

1 

-1 

1 

0.1 

0.2 

1 

3.9 

3.9 

500 

1 

E 

6.9 

lJJ 

-1 

O.Il 

1 

3.9 

3.9 

500 

I 

F’igure  6.5  displays  the  constant  voltage  contours  using  the  data  for  case  A.  The  left  electrode 
is  at  a  1  volt  potential,  while  the  right  electrode  and  bottom  border  are  grounded  {V  =  0).  All 


6-7 


Constant  Voltage  Contours 


-1.000  -0.600  -0.200  0.200  0.600  1.000 

X  position 


Figure  6.5.  Constant  voltage  contours  in  cell  using  the  case  A  data  set 

layers  contain  equal  permittivity  (cr  =  !)•  The  curves  are  smooth  across  both  boundaries  (layers 
2-3  and  layers  3-4),  which  demonstrates  that  the  numerical  approximation  of  the  solution  matches 
very  well  at  the  boundaries.  The  plot  contains  100  voltage  levels,  so  each  curve  represents  a  step 
of  0.01  volts. 

Figure  6.6  shows  the  constant  voltage  contours  using  the  data  for  case  B.  The  left  and  right 
electrodes  have  1  and  -1  volt  applied  respectively,  while  the  bottom  border  is  grounded.  All  layers 
contain  equal  permittivity  (e^  =  !)■  Again  the  curves  are  smooth  across  the  boundaries  and  100 
contour  levels  are  shown,  but  the  levels  are  now  incremental  steps  of  0.02  volts  (2  volt  difference 
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Figure  6.6.  Constant  voltage  contours  in  cell  using  the  case  B  data  set 

between  electrodes).  The  strong  electric  field  near  the  electrode  corners  is  clearly  illustrated,  and 
the  even  symmetry  is  also  apparent  in  the  plot. 

The  constant  voltage  contours  using  the  data  for  case  C  is  shown  in  Figure  6.7.  The  electrodes 
are  set  to  1  and  -  I  volts,  with  the  bottom  border  grounded.  However,  a  high  permittivity  layer 
(fr  =  100)  has  been  introduced  between  the  electrodes.  The  layer  is  well  defined  and  the  strong 
electric  field  within  the  layer  is  observed  in  the  contour  plot.  Again,  100  contour  levels  are  shown 
with  steps  of  0.02  volts  between  curves. 

Case  D  (Figure  6.8)  introduces  a  thick,  high  permittivity,  layer  4  (fr  =  500),  while  layers  2 
and  3  are  insulating  material  (e^  =  3.9)  and  layer  5  is  air  (fr  =  1).  The  electrodes  have  applied 


U  position 


Figure  6.8.  Constant  voltage  contours  in  cell  using  me  case  D  data  set 


potentials  of  1  and  -1  volt,  with  the  bottom  of  layer  2  grounded.  Layer  4  is  well  defined  in  the  plot 
and  a  strong  electric  field  is  observed  within  the  layer.  Again  100  contour  levels  are  displayed  with 
a  step  of  0.02  volts  between  curves.  The  solution  is  smooth  across  layers  2  and  3  since  they  contain 
equal  permittivities. 

Figure  6.9  and  6.10  display  the  constant  voltage  contours  of  the  sensor  cell  with  parameters 
selected  to  represent  the  CHEMFET.  Layer  2  is  a  thin  silicon  dioxide  layer  (cr  =  3.9),  layer  3  is 
F^-glass  (fr  =  3.9),  layer  4  is  a  very  thin  polymer  layer  (fr  =  500),  while  layer  5  is  air  (fr  =  1).  The 
bottom  border  is  grounded  and  the  electrodes  are  at  ±1  volt.  Layer  4  is  not  easily  seen  due  to  the 


Constant  Voltage  Contours 


S  xlO° 


Figure  6.9.  Constant  voltage  contours  of  the  sensor  ceU  with  *10  contour  levels 

scale  of  the  plot.  Figure  6.9  contains  *10  contour  levels  with  a  step  of  0.05  volts  between  curves, 
while  Figure  6.10  contains  80  contour  levels  with  steps  of  0.025  volts  between  curves. 

The  contour  plots  shown  in  Figures  6.5-6.10  illustrate  that  the  solution  behaves  as  expected 
and  provide  a  visual  representation  of  the  solution  under  different  conditions.  Boundaries  between 
layers  are  unrecognizable  if  the  permittivities  are  identical  across  the  interface.  When  the  permit¬ 
tivities  do  not  match,  the  boundary  becomes  welt  defined. 

6.1.3  Capacitance  Graphs.  With  approximations  for  the  potential  coefficients,  the  charge, 
and  the  capacitive  elements  can  be  calculated,  using  Equations  -1.30,  4.60,  and  4.63.  Thus,  we 
can  experiment  with  the  results  and  observe  the  change  in  the  capacitive  elements  if  we  vary 
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key  parameters  in  the  periodic  cell.  This  investigation  will  examine  the  interdigitated  electrode 
structure  with  and  without  the  applied  sensor  coating.  The  model  parameters  established  as 
the  baseline  for  the  investigation  are  shown  in  Table  6.4.  These  values  were  selected  using  the 
CHEMFET  structure  as  a  guide  and  the  cell  dimensions  were  normalized  to  the  electrode's  width. 
The  cell  capacitance  terms  are  all  expressed  on  a  per  unit  length  basis;  that  is,  with  units  of  Farads 
per  centimeter  (F/cm). 

Table  6.4.  Normalized  model  parameters  used  as  the  baseline  for  capacitance  calculations 


material  above  electrodes 

wm 

mm 

1 

D 

B9 

in 

n 

ESI 

air 

0.11 

3 

3.9 

1 

1 

ri 

insulator/sensor  coating 

0.1 

0.1 

0.11 

Ji 

3.9 

3.9 

500 

Bl 

6. 1.3.1  Effects  of  Insulator  Thickness  on  the  Capacitive  Elements.  Figure  6.11  dis¬ 
plays  the  capacitive  elements,  Cg  and  Ce,  relative  to  a  range  of  insulator  thicknesses,  with  the 
remaining  parameters  set  to  their  baseline  values.  The  capacitance  between  electrodes,  Ce,  in¬ 
creases  with  insulator  thickness,  while  the  capacitance  to  ground,  Cj,  decreases.  Figure  6.12  shows 
similar  results  for  the  case  when  the  sensor  coating  is  considered.  Again,  Cg  decreases  with  increas¬ 
ing  insulator  thickness,  while  Ce  increases.  While  Cg  increases  slightly  when  the  sensor  coating  is 
considered,  Ce  increeises  by  a  factor  of  10.  To  understand  these  results,  an  analysis  of  the  cell  is 
needed. 

With  a  thin  insulating  layer,  most  of  the  electric  flux  exists  between  each  electrode  and  the 
ground  plane.  As  the  insulator  thickness  is  increased,  there  will  be  a  proportionately  greater  flux 
between  the  electrodes  and  a  lesser  amount  between  each  electrode  and  the  ground  plane.  Therefore, 
as  the  insulator  thickness  is  increased,  we  should  expect  Ce  to  increase  and  Cg  to  decrease. 


6. 1.3. 2  Effects  of  Electrode  Height  on  the  Capacitive  Elements.  Figure  6. 1 3  shows  the 
capacitive  elements,  Cg  and  C'e,  relative  to  a  range  of  electrode  heights,  with  the  remaining  model 
parameters  set  at  their  baseline  values.  In  the  case  of  air  above  the  electrodes,  the  numerical  results 


Figlire  6.1 1.  Capacitive  element  graphs  for  varying  d,  with  air  above  the  electrodes 
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(a)  Cell  insulator  capacitance  per  unit  electrode  length 
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(b)  Cell  inter-terminal  capacitance  per  unit  electrode  length 
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Figure  6.12.  Capacitive  element  graphs  for  varying  d,  with  the  sensor  coating  considered 
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shows  that  both  Cg  and  Ce  increase  with  electrode  height.  Figure  6.14  displays  similar  results  for 
the  capacitive  elements  with  the  sensor  coating  considered.  Again,  Ce  increases  with  electrode 
height,  but  now  Cg  decreases.  Also,  Ce  increases  by  a  factor  of  10  when  the  sensor  coating  is 
considered.  To  understand  the  influence  of  the  electrode  height  on  the  capacitive  elements,  an 
analysis  of  the  cell  is  required. 

Increasing  the  electrode  height  provides  more  surface  area  for  charge  storage,  which  means 
that  the  electric  field  between  the  electrodes  increases.  Therefore,  increasing  the  electrode  height 
should  increase  Ce-  With  air  above  the  electrodes,  the  potential  functions  across  the  layer  3  inter¬ 
faces  are  non-linear  (see  Figure  6.3),  indicating  that  most  of  the  cell’s  electric  flux  is  concentrated 
between  the  electrodes  and  the  ground  plane  rather  than  between  the  two  electrodes.  For  this  case, 
increasing  the  electrode  height  should  increase  Cg  due  to  the  increased  charge  storage  capacity. 
With  the  sensor  coating  considered,  the  potential  functions  along  the  layer  3  interfaces  are  sig¬ 
nificantly  more  linear  (see  Figure  6.4),  indicating  that  more  flux  is  concentrated  between  the  two 
electrodes.  In  this  case,  increasing  the  electrode’s  height  will  decrease  Cg  due  to  the  correspondingly 
increased  electric  field  between  the  two  electrodes. 

6. 1.3. 3  Effects  of  Electrode  Separation  on  the  Capacitive  Elements.  Figure  6.15  dis¬ 
plays  the  capacitive  elements,  Cg  and  Ce,  relative  to  a  range  of  electrode  separations,  with  the 
remaining  model  parameters  set  to  their  baseline  values.  The  numerical  results  show  that  Cg  in¬ 
creases  with  larger  separation,  approaching  the  limiting  value  for  the  parallel  combination  of  two 
tincoupled  micro  .ips  in  the  layered  materials.  The  results  also  indicate  that  the  inter-terminal 
capacitance,  Ce,  decreases  with  larger  separation.  Figure  6.16  shows  similar  results  for  the  ca¬ 
pacitive  elements  when  the  sensor  coating  is  considered.  Again,  Cg  increases  with  larger  electrode 
separation  to  a  limiting  value,  while  Ce  decreases.  To  understand  these  results,  an  analysis  of  the 
cell  is  needed. 
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(a)  Cell  insulator  capacitance  per  unit  electrode  length 
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(b)  Cell  inter-terminal  capacitance  per  unit  electrode  length 
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Figure  6.13.  Capacitive  element  graphs  for  varying  /i,  with  ai'  above  the  electrodes 
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(a)  Cell  insulator  capacitance  per  unit  electrode  length 
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Figure  6.14.  Capacitive  element  graphs  for  varying  h,  with  the  sensor  coating  considered 
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Increasing  the  electrode  separation  decreases  the  electric  field  between  the  electrodes  and  this 
phenomenon  magnifies  the  effect  of  the  ground  plane  on  each  conductor.  Therefore,  the  capaci¬ 
tance  to  ground  should  increase  with  more  separation  until  the  cell  is,  effectively,  two  uncoupled 
microstrips,  while  the  capacitance  between  the  electrodes  should  de  -ase  with  more  separation. 
In  addition,  the  extreme  case  of  s  =  0  will  result  in  the  cell  effectively  having  only  one  electrode, 
with  the  capacitance  to  ground  becoming,  essentially,  a  parallel  plate  capacitance. 

6. 1.3.4  Effects  of  Sensor  Coating  Thickness  on  the  Capacitive  Elements.  Figure  6.17 
shows  the  capacitive  elements,  Cg  and  Ce,  relative  to  a  range  of  sensor  coating  thicknessess,  with  the 
remaining  cell  parameters  set  to  their  baseline  values.  The  numerical  results  show  that  an  increase 
in  sensor  coating  thickness,  p,  will  increase  both  Cg  and  Ce-  From  the  graphs,  Cg  is  more  sensitive 
to  changes  in  the  sensor  coating  thickness  compared  to  Ce  when  the  coating  thickness  is  at  the 
thinner  end  of  the  range.  However,  Ce  is  more  sensitive  to  changes  in  the  sensor  coating’s  thickness 
compared  to  Cg  when  the  coating  thickness  is  at  the  thicker  end  of  the  range.  To  understand  the 
overall  effect  of  the  sensor  coating  thickness  on  the  capacitance  values,  an  analysis  of  the  cell  is 
necessary. 

The  sensor  coating  layer  has  a  large  relative  permittivity,  and  its  presence  was  observed  to 
affect  the  electic  flux  of  the  periodic  cell  (see  Figures  6.3  and  6.4).  The  presence  of  the  sensor  coating 
acts  to  increase  the  electric  field  strength  between  the  electrodes,  and  it  correspondingly  reduces 
the  influence  of  the  ground  plane  on  the  conductors.  Since  the  insulator  is  so  thin  (the  electrodes 
are  much  closer  to  the  ground  plane  compared  to  their  separation),  some  of  the  electric  flux  within 
the  sensor  coating  reaches  the  ground  plane.  Therefore,  increasing  the  sensor  coating  thickness 
should  increase  both  the  capacitance  between  the  electrodes  and  the  capacitance  to  ground. 

6. 1.3. 5  Effects  of  Sensor  Coating  Permittivity  on  the  Capacitive  Elements.  Figure  6.18 
displays  the  numerical  results  of  the  capacitive  elements,  Cg  and  Ce,  relative  to  a  range  of  relative 
permittivities  for  the  sensor  coating,  with  the  remaining  model  parameters  set  to  their  baseline  val- 
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Figure  6.15.  Capacitive  element  graphs  for  varying  s,  with  air  above  the  electrodes 


Figure  6.17,  Capacitive  element  graphs  for  varying  p,  with  the  sensor  coating  considered 


lies.  The  results  show  that  both  Cg  and  Ce  increase  with  an  increase  in  the  relative  permittivity 
of  the  sensor  coating  (£4  =  where  f„  is  the  permittivity  of  free  space).  These  results  Indicate 

that  Cg  is  strongly  dependent  on  the  sensor  coating’s  1  laterial  properties  (this  was  not  expected  at 
the  outset  of  this  investigation)  and  provide  evidence  that  both  capacitive  elements  play  a  major 
role  in  the  sensing  action  of  such  a  device.  As  discussed  in  the  previous  section,  the  sensor  coating 
layer  affects  both  Cg  and  Ce-  Therefore,  an  increase  in  the  sensor  coating’s  permittivity  should 
increase  both  'apacitive  elements. 

6.2  Results  of  the  Resistive  Solution 

A  FORTRAN  program  was  created  to  solve  the  matrix  equation  in  section  5.4,  and  thereby 
approximate  the  potential  within  the  resistive  layer.  The  program  utilized  the  IMSL  gaussian 
elimination  routine  “LSARG”  [7:11].  The  results  were  then  used  to  calculate  the  resistance  of  the 
sensor  coating  layer  with  varying  parametric  conditions,  using  Equation  5.22. 

6.2.1  Ulustrating  the  Resistive  Potential  Solution.  To  illustrate  the  solution,  the  potential 
along  the  lower  boundary  is  displayed  in  Figure  6.19.  The  displayed  potential  satisfies  the  boundary 
data  of  the  problem;  specifically,  the  potential  on  the  left  electrode  is  +0.5  volts  while  the  potential 
on  the  right  electrode  is  —0.5  volts.  The  potential  is  essentially  linear  between  the  electrodes.  To 
illustrate  the  potential  distribution  within  the  layer,  it  v  as  calculated  at  discrete  points  on  a  100 
by  50  point  grid.  The  grid  spacing  in  the  x  direction  is  (s  +  1)/100,  while  the  spacing  in  the  y 
direction  is  p/50.  A  plot  of  constant  voltage  contours  was  created  using  AFIT’s  METALIB  plotting 
routines  [1].  Figure  6.20  displays  the  constant  voltage  contour  plot  (50  contour  levels)  with  steps  of 
0.02  volts  between  curves.  Figure  6.21  shows  the  constant  voltage  contours  of  the  same  grid  with 
100  contours  and  steps  of  0.01  volt  between  curves. 

6.2.2  Resistive  Element  Graphs.  With  estimates  for  the  coefficients  of  the  potential  func¬ 
tion.  the  resistance  of  the  layer  was  determine  1  using  Equation  5.22.  The  behavior  of  the  resistive 
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element  can  be  investigated  by  experimenting  with  the  solutions  and  results  of  the  resistive  solution 
while  varying  key  parameters  of  the  periodic  cell.  All  resistance  values  displayed  are  expressed  on 
a  unit  length  basis,  with  units  of  Ohms  centimeter  (O-cm).  The  baseline  model  parameters  used 
to  investigate  the  resistive  element  are  shown  in  Table  6.5.  Again,  these  cell  parameter  values  were 
selected  using  the  CHEMFET  structure  as  a  guide  and  the  cell  dimensions  were  normalized  with 
respect  to  the  electrode’s  width. 

Table  6.5.  Normalized  model  parameters  used  as  baseline  for  the  resistive  calculations 


parameter 

P 

(T4 

1  value 

0.01 

Tj 

10-® 

6.2.2. 1  Effects  of  Sensor  Coating  Thickness  on  the  Resistive  Element.  Figure  6.22 
shows  the  behavior  of  the  model  resistive  element.  /Ip,  relative  to  a  range  of  values  for  the  sen.sor 
coating  thickness,  with  the  remaining  parameters  set  to  their  baseline  values.  These  results  show 
that  as  the  sensor  coating’s  thickness  increases,  the  resistance  decreases.  With  a  thicker  sensor 
coating  layer,  there  is  more  resistive  material  present;  and  hence,  the  resistance  through  the  material 
should  decrease.  In  the  CHEMFET  application,  the  sensor  coating  is  expected  to  be  thin  (small 
relative  to  the  other  cell  dimensions,  including  the  oxide  thickness),  since  the  coating  is  applied  by 
a  vapor  deposition  process. 

6. 2. 2. 2  Effects  of  Electrode  Separation  on  the  Resistive  Element.  Figure  6.23  displays 
the  behavior  of  the  resistive  model  element,  Rp,  relative  to  a  range  of  the  electrode's  separation, 
with  the  remaining  model  parameters  set  to  their  baseline  values.  With  a  thin  sensor  coating  layer 
(p  =  0  01),  there  is  a  linear  relationship  between  the  electrode  separation  and  Rp.  This  linear 
r<'l.'itionship  permits  a  coni[iari.son  with  the  results  obtained  from  the  formula  for  a  parallel  plate 
ri.si.stor  approximation  Witli  this  a-ssumption,  the  inter-terminal  resistance  for  a  imit  length.  Rp. 
Ix'f  ()rnes 
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Figure  6.22.  Cell  resistance  graph  for  varying  p 


Rp  = 


s 

(TP 


(G.l) 


Using  the  baseline  values,  the  parallel  plate  approximation  yields  1  Gfi-cm,  while  the  numerical 
solution  is  1.01  GQ-cm.  Therefore,  the  parallel  plate  approximation  matches  the  numerical  solution 
to  within  1%  for  this  one  set  of  calculations,  which  considered  a  thin  layer  (p  <C  .s). 


6. 2. 2. 3  Effects  of  Sensor  Coating  Conductivity  on  the  Resistive  Element.  Figure  6.2'4 
shows  the  effect  of  the  sensor  coating’s  conductivity  on  the  model  element  Rp.  These  results  demon¬ 
strate  a  strong  dependence  of  the  resistance  element  to  the  sensor  coating's  material  properties,  as 
expected  from  the  parallel  plate  estimate.  Therefore,  any  change  in  the  sensor  coating’s  cori'*"c- 
tivity  directly  affects  the  resistance  between  the  electrodes,  and  suggests  that  the  resit,  ..lent 

is  also  a  major  factor  in  the  sensing  .action  of  the  structure. 


s,  normalized  to  the  electrode’s  width 


Figure  6.23.  Cell  resistance  graph  for  varying  s 
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Figure  6.24.  Cell  resistance  graph  for  varying  conductivity  in  the  sensor  coating 


VII.  Frequency  Dependence  of  the  .Model  Elements 


In  this  chapter,  an  analysis  of  the  complex  boundary  condition  at  the  interfaces  between 
the  material  layers  is  used  to  determine  the  frequency  dependence  of  the  model  elements.  Thus 
far,  the  analysis  leading  to  the  capacitive  and  resistive  elements  was  based  upon  the  high  and 
low  frequency  arguments  of  Equation  3.41.  Therefore,  the  previous  results  were  limited  to  the 
high  frequency  capacitance  and  DC  (direct  current)  resistance.  It  would  be  beneficial  to  verify  the 
validity  of  the  previous  results  and  determine  the  frequency  dependence  of  the  model  elements. 
Towards  this  objective,  let  us  consider  the  complex  boundary  condition  presented  in  Equation  3.41 

(7.1) 

where  (  is  the  permittivity  of  the  material,  a  is  the  conductivity  of  the  material,  w  is  the  angular  ex¬ 
citation  frequency,  and  Eni.  is  the  normal  conponent  of  the  electric  field  in  layer  k.  dimensionless 
complex  relative  permittivity,  e',  can  be  defined  as 

(7.2) 

WCo 

r 

where  is  the  permittivity  of  free  space  and  tr  is  the  real  relative  permittivity  of  the  material 
(f  =  frfo).  Using  Equation  7.2  in  7.1  yields 

'■1  (7 .3) 

Using  Equation  7.3,  the  capacitive  (or  high  frequency)  problem  formulation  can  be  extended  to  all 
frequencies  by  replacing  the  real  relative  permittivity  with  the  complex  formulation,  provided  the 
potential  is  also  complex  so  as  to  satisfy  the  complex  boundary  conditions. 

Using  the  complex  relative  permittivity  defined  in  Equation  7.2  and  the  complex  I.MSL  gaus- 
sian  elimination  routine  “LSACG”  [7:31],  the  complex  coefficients  of  the  potential  in  each  region 
can  be  found.  With  the  potential  known,  the  electric  field  can  be  determined  which  leads  to  both 
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the  real  charge  on  the  electrode,  as  well  as  the  real  current  flowing  between  the  terminals.  From 
the  charge  and  current,  the  capacitance  is  found  using  Equations  4.60  and  4.63.  The  resistance  is 
found  using  Equations  4.3,  5.14,  5.16  and  integrating  the  resulting  current  density,  as  implied  in 
Equation  5.20,  which  yields 

Rp  =  (743?  I^dnsin  exp(-A4/)  ~  j  | 

where  3?  denotes  the  real  part  and 


—  sinh[A4(/ — /i)  +  —  cosh(A4(/ — /i)]  ('-6) 

Therefore,  the  model  elements  can  be  determined  for  any  given  frequency  of  operation. 

Figure  7.1  displays  the  model  elements  (Cj,  Ce,  and  Rp)  as  a  function  of  frequency  using 
the  model  data  shown  in  Table  7.1.  As  in  the  previous  chapter,  the  cell  parameters  were  selected 
using  the  CHEMFET  structure  as  a  guide,  and  the  cell  dimensions  were  normalized  with  respect 
to  the  electrode’s  width.  In  addition,  the  model  elements  are  all  expressed  on  a  unit  length  basis, 
with  the  capacitive  elements  values  in  Farads  per  centimeter  (F/cm)  and  resistance  values  in  ohms- 
centimeter  (Q-cm). 


Table  7.1.  Normalized  data  used  for  the  frequency  analysis 


parameter 

d 

wm 

1 

m 

BH 

B9 

n 

EH 

0-2 

(74 

0-5 

value 

0.1 

0.1 

I  0.11 

3.9 

3.9 

500 

1 

0 

■H 

ISBI 

Bli 

Throughout  the  range  of  frequencies  shown,  the  sensor  coating  behaves  both  as  a  good  con¬ 
ductor  and  as  a  lossy  dielectric.  The  transition  between  these  two  modes  is  illustrated  in  Figure  7.1, 
atul  it  occurs  when 
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Figure  7.1.  Behavior  of  the  model  elements  versus  frequency 


where  Wp  is  the  angular  transition  frequency  of  the  sensor  coating,  which  was  discussed  in  Chap¬ 
ter  III  and  defined  as  the  frequency  where  the  sensor  coating  transitions  from  a  conductor  to  a 
lossy  dielectric.  From  the  data  in  Table  7.1,  the  angular  transition  frequency  is 


cta  10“® 

ujp  =  —  =  - =  226.2  rad/sec 

^  u  SOOfo  ' 

which  leads  to  an  ordinary  transition  frequency  (/p)  of 
/p  =  ^r.36  Hz 

aTT 


(7.8) 


(7.9) 


This  result  is  clearly  in  agreement  with  Figure  7.1.  Element  values  far  below  the  transition  are 
suspect.  The  matrix  which  leads  to  the  potential  coefficients  has  a  large  condition  number  (10®) 
at  the  lowest  frequency  shown,  and  the  IMSL  routine  warns  of  ill-conditioning  which  suggests  that 
the  result  may  be  inaccurate.  This  warning  is  given  for  any  frequency  of  operation  below  3  Hz  with 
the  data  shown  in  Table  7.1.  The  ill-conditioning  flag  is  dependent  on  the  numerical  precision  of 
the  computer.  The  FORTRAN  programs  were  processed  on  an  ELXSI  System  6400  using  single 
precision.  The  program  was  also  processed  using  double  precision,  which  resulted  in  identical 
output  but  no  warnings. 

In  addition,  the  low  frequency  resistance  value  shown  is  approximately  10.6  Gfi-cm,  while  the 
calculated  value  using  the  resistive  solution  is  10.1  Gfi-cm.  Figure  7.2  displays  resistance  versus 
frequency  behavior  using  50,  70,  100,  150,  and  200  coefficients  for  the  potential.  Since  the  potential 
solution  is  a  convergent  infinite  series,  more  coefficients  should  yield  a  better  approximation  to  the 
solution;  however,  the  processing  lime  increases  significantly  to  obtain  these  additional  coefficients 
(for  example  it  takes  15.7  CPU  seconds  to  determine  50  coefficients,  25.3  sec  for  70,  64.1  sec  for 
100,  194.3  sec  for  150,  and  423.9  .sec  for  200).  With  more  coefficients,  the  low  frequency  resistance 
becomes  smaller;  with  200  coefficients,  the  low  frequency  resistance  is  approximately  10.2  Gfl-cm. 
However  the  ill-conditioning  warning  occurs  at  slightly  higher  frequencies  when  more  terms  are 
u.sed,  and  with  200  coefficients,  the  IMSL  routine  provides  warnings  for  frequencies  below  15  Hz 


Figure  7.2,  Rp  versus  frequency  (50  through  200  coefficients) 


(below  4  Hz  for  70,  8  Hz  for  100,  and  11  Hz  for  150).  Therefore,  the  data  is  suspect  for  frequencies 
less  than  fp.  Figures  7.3  and  7.4  provide  similar  data  for  the  model  elements  Cg  and  C^. 


7.t  Effects  of  Conductivity  on  the  Frequency-Dependent  Model  Elements 

Figure  7.5  displays  the  model  element  Cg  versus  frequency  for  —  10“',  10“^,  and  10“^° 
15/cm,  with  the  other  model  parameters  shown  in  Table  7.1.  Figures  7.6  and  7.7  depict  similar 
results  for  the  model  elements  Ce  and  ftp.  For  each  model  element,  the  curves  are  identical  but 
shifted  in  frequency,  which  clearly  illustrates  tne  elements’  dependence  on  the  sensor  coating’s 
conductivity.  For  each  conductivity  value,  an  ill-conditioned  warning  was  returned  for  the  lowest 
frequency  shown,  which  is  consistent  with  the  prior  results  and  suggests  a  numerical  instability  at 
the  lower  frequencies. 
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1.2  Effects  of  Permittiinty  on  the  Frequency-Dependent  Model  Elements 


Figure  7.8  displays  for  =  20,  500,  and  1300,  with  the  other  model  data  shown  in 
Table  7.1.  Figures  7.9  and  7.10  depict  similar  graphs  for  Cf;  and /?p.  While  the  curves  for  tr,  =500 
and  1300  appear  consistent  with  the  results  for  the  various  conductivities  considered,  this  is  not 
the  case  for  £r<  =  20.  For  Cj,  with  =  20,  the  curves  again  appear  i(.entical  but  shifted;  however, 
this  illusion  is  incorrect.  For  (r,  =  20,  the  sensor  coating  transition  frequency  (fp)  is  900  Hz,  and 
the  curve  should  begin  to  depart  from  the  high  frequency  value  at  approximately  6000  Hz  (see  C'e 
and  Rp  graphs  for  comparison);  however,  it  appears  to  deviate  from  the  high  frequency  value  at 
approximately  300  Hz,  which  is  where  Rp  reaches  an  intermediate  peak  between  the  high  and  low 
frequency  values.  In  addition,  Cg  becomes  negative  (inductive?)  at  this  same  point.  For  Cr,  =  20, 
no  ill-conditioned  warnings  were  received  for  the  frequencies  shown,  while  ill-conditioning  warnings 
were  given  for  the  first  frequency  shown  for  €r«  =  500  and  1300. 

For  =  20,  frequencies  below  fp  do  not  behave  as  expected  for  both  Ce  and  Rp.  While  the 
general  shape  of  Ce  versus  frequency  appears  valid  (for  frequencies  above  100  Hz),  the  capacitance  is 
not  expected  to  become  negative.  This  could  be  an  indication  of  inductive  behavior,  or  it  could  also 
be  another  indication  that  the  element  values  predicted  for  lower  frequencies  are  suspect  (possibly 
invalid).  For  Rp,  the  intermediate  peak  between  the  low  and  high  frequency  values  is  interesting. 
It  could  signify  a  possible  resonance,  or  again  indicate  the  data  is  suspect  at  the  lower  frequencies. 
For  comparison,  graphs  of  the  element  values  versus  frequency  are  shown  for  fp,  =  IOOt  75,  and  50 
in  Figures  7.11,  7.12,  and  7.13.  For  fp,  =  100,  an  intermediate  peak  for  Rp  is  observed,  and  Ce 
becomes  negative  for  the  lower  frequencies.  As  the  relative  permittivity  is  reduced,  the  Rp  peak 
becomes  more  pronounced,  and  Ce  becomes  more  negative  at  the  lower  frequencies. 

For  frequencies  above  fp,  there  is  good  agreement  with  all  the  element  graphs,  and  the  high 
frequency  capacitance  values  match  exactly  the  capacitance  values  found  using  the  capacitive  (high 
frequency)  solution.  The  high  frequency  resistance  can  be  determined  by  treating  the  sensor  coating 
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as  a  lossy  dielectric  with  the  limit  taken  as  w  —  oo.  The  current  through  the  sensor  coating  layer 
can  be  found  by  integrating  the  electric  field  of  the  capacitive  solution  as  was  described  for  the 
resistance  calculation  in  Chapter  V.  The  resistance  is  then  calculated  from  Ohm’s  law. 

1.3  Comparison  of  the  Model  Element  Values  (High  Frequency  and  Frequency- Dependent) 

The  graphs  of  the  capacitive  model  elements  versus  frequency  illustrated  that  the  element 
values  attained  a  high  frequency  limit;  these  values  agree  with  the  capacitive  values  found  using  the 
high  frequency  boundary  conditions.  The  high  frequency  resistance  was  determined  by  considering 
the  sensor  coating  as  a  lossy  dielectric  in  the  limit  as  w  — •  oo  (where  =  f^),  and  calculating 
the  current  through  the  layer.  While  the  complex  permittivity  solution  leads  to  the  more  correct 
model  element  values  for  any  finite  frequency,  this  approach  also  increases  the  computational  time 
needed  to  find  the  element  values.  Hence,  if  we  use  the  high  frequency  values,  we  can  reduce  the 
computation  time  and  still  have  reasonably  accurate  lumped  element  values  for  a  given  frequency 
range.  To  compare  the  high  frequency  values  relative  to  the  frequency-dependent  values,  the  input 
and  inter- terminal  impedance  are  considered. 

The  inter-terminal  impedance,  Zt,  is  the  parallel  combination  of  Rp  and  Ce-  Since  admit¬ 
tances  in  parallel  add,  the  inter-terminal  admittance,  Yt,  is 

YT=Gp-^juCE  (7.10) 

where  Yt  =  1/^t,  Gp  =  If  Rp,j  =  1,  and  a;  is  the  angular  excitation  frequency.  The  admittance 

to  ground,  Yg  is 

yS=j^Cg  (7.11) 

Let  Zfj  denote  the  impedance  of  the  series  combination  of  Zt  and  Zg,  where  Zg  =  l/Yg.  Then 

ZN  =  ZT+Zg  (7.12) 
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(a)  frequency-dependent  Ce,  with  fr4  =  20 
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(b)  frequency-dependent  Ce,  with  (ri  =  500  and  (T4  =  10“® 
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(c)  frequency-dependent  Ce,  with  Cr4  =  1300 


Figure  7.9.  Ce  versus  frequency  for  Cr,  =  20,  500,  1300 
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(a)  frequency-dependent  Cg,  with  Cr4  =  100 


(b)  frequency-dependent  Ce,  with  CrA  =  100 


(c)  frequency-dependent  Rj,,  with  er4  =  100 


Figure  7.11.  Model  Elements  versus  frequency  for  tr,  =  100 
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Figure  7.12.  Model  Elements  versus  frequency  for  fr,  =  75 


The  input  impedance,  Z,n,  is  liie  impedance  at  the  input  terminals  of  the  device,  which  is  the 
parallel  combination  of  and  Zg.  Using  the  fact  that  admittances  in  parallel  add,  yields 

v;„  =  y;v  +  v;,  (7.13) 

where  V',„  =  XfZin- 

Figure  7.14  displays  the  input  impedance  of  the  interdigitated  electrode  sensor  structure 
(both  magnitude  and  phase),  using  both  the  frequency-dependent  model  element  values  and  the 
high  frequency  values.  The  two  curves  are  identical  at  the  higher  frequencies,  and  they  only  deviate 
at  the  lower  frequencies  (as  expected).  Even  at  the  lower  frequencies  (near  /p),  the  input  impedance 
has  only  a  small  deviation  (in  botii  magnitude  and  phase)  relative  to  using  the  different  values  for 
the  model  elements. 

Figure  7.15  displays  the  inter-terminal  impedance  of  the  interdigitated  electrode  sensor  struc¬ 
ture  (both  in  magnitude  and  phase),  relative  to  using  both  the  frequency-dependent  values  and 
the  high  fretjuency  values.  Again  the  two  curves  are  nearly  identical,  with  only  a  relatively  small 
deviaticm  at  the  lower  frequencies. 

In  conclusion,  the  model  element  values  found  using  the  high  frequency  (capacitive)  boundary 
conditions  provide  a  valid  model  for  /  >  fp.  While  more  accurate  model  element  values  can  be 
obtained  using  frequency  analysis  with  a  complex  dielectric  constant,  tiie  additional  computaticn 
time  may  not  be  beneficial. 


7-18 


Figure  7.14.  Input  impedance  using  frequency-dependent  and  high  frequency  elements  values 


(a)  inter-terminal  impedance  for  a  unit  length  -  magnitude 
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(b)  inter-terminal  impedance  for  a  unit  length  -  phase 
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VIH.  Conclusion 


S.l  Summary 

A  lumped-parameter  electrical  model  for  the  interdigitated  electrode  sensor  structure  has 
been  developed.  Three  elements  have  been  identified;  they  are  a  capacitance  from  each  of  the  two 
electrodes  relative  to  ground  (Cg),  a  capacitance  between  electrodes  (Cfr),  and  a  resistance  through 
the  sensor  i-oating  (Rp).  The  elements  were  determined  from  the  electric  field  within  the  structure 
by  first  solving  for  the  potential  distribution. 

•A  periodic  cell  was  established  to  represent  the  configuration  of  electrodes,  since  the  inter¬ 
digitated  electrode  sensor  structure  contains  a  repetitive  pattern  in  the  plane  transverse  to  the 
structure’s  length.  Five  layers  within  the  cell  were  then  defined.  The  interdigitated  electrode  .sen¬ 
sor  structure  was  then  reduced  to  a  two-dimensional  problem  using  a  transverse  electromagnetic 
(TE.M)  wave  analysis  and  limiting  the  model  to  a  frequency  range  where  the  potential  is  e.ssen- 
tially  constant  along  the  structure’s  length.  This  situation  also  resulted  in  the  potential  satisfying 
baplace’s  equation  in  the  transverse  plane.  Boundary  conditions  at  ‘he  interfaces  between  layers 
were  then  developed  for  the  low  and  high  frequency  cases,  which  correspond  to  the  conditions  that 
the  sensor  coating  is  either  a  good  conductor,  or  a  lossy  dielectric,  respectively. 

The  capacitive  elements  were  determined  using  boundary  conditions  on  the  potential  based 
on  the  high  frequency  limit.  A  Fourier  series  representation  of  the  potential  weis  established  for  each 
layer  in  the  cell,  with  the  Fourier  coefficients  determined  by  the  boundary  conditions.  These  coeffi¬ 
cients  were  approximated  numerically  by  truncating  the  infinite  Fourier  series.  From  the  potential, 
the  electric  field  was  found,  along  with  the  charge  enclosed  on  an  electrode.  The  capacitance  values 
were  determined  from  the  charge  calculations  using  Weeks’  method  [17]. 

The  resistive  element  was  determined  using  boundary  conditions  for  the  potential  based  on  a 
low  frequency  limit.  Again,  a  Fourier  .series  representation  of  the  potential  was  established  for  each 
region  within  the  layer,  and  the  Fourier  coefficients  were  obtained  from  the  boundary  conditions 


at  the  interfaces  between  regions.  These  coefficients  were  also  estimated  by  truncating  the  infinite 
Fourier  series.  From  the  potential,  the  electric  field  and  current  density  were  found  through  the 
partially  tonductive  sensor  coating  layer.  The  resistance  was  finally  calculated  using  Ohm’s  law 

The  capacitive  and  resistive  solutions  were  illustrated  using  specific  cell  parameter  values  that 
would  highlight  certain  aspects  of  the  solution.  At  this  point,  the  investigation  shifted  from  a  general 
nature  to  a  particular  case  focusing  on  i.he  dimensional  and  material  values  closely  related  to  the 
chemically-sensitive  field-effect  transitor  (CHEMFET1.  Then,  the  capacitive  and  resistive  elements 
were  calculated  for  a  range  of  structural  dimensions  sense  r  coating  permittivities,  and  sensor 
coating  conductivities.  The  results  were  pr<"sented,  and  they  revealed  that  all  the  model  elements 
were  sensitive  to  the  variable  electrical  properties  of  the  thin  sensor  coating  layer  covering  the 
electrodes.  In  addition  the  dependence  of  the  model  elements  on  each  of  the  structural  dimensions 
and  material  properties  was  illustrated. 

Using  a  complex  permittivity,  the  frequency-dependent  lumped  parameters  were  obtained; 
these  r‘=sults  agreed  with  the  high  frequency  capacitive,  and  low  frequency  resistive  values  in  the 
appropriate  limits.  From  the  frequency  plots  shown  in  Chapter  VII,  the  transition  from  the  high  fre¬ 
quency  to  the  low  frequency  values  was  illustrated,  corresponding  to  the  sensor  coating’s  transition 
from  a  conductor  to  a  lossy  dielectric.  For  frequencies  larger  than  the  sensor  coating’s  transition 
point,  the  numerical  results  showed  that  all  element  values  approach  a  high  frequency  value  over  a 
range  of  sensor  coating  permittivties  and  conductivities.  At  frequencies  below  the  transition,  the 
numerical  results  are  suspect  due  to  a  possible  numerical  instability  in  the  solution.  In  general, 
the  results  at  the  lower  frequencies  do  not  behave  as  expected.  The  potential  solution  for  the  DC 
(direct  current)  case  must  satisfy  different  boundary  conditions  at  the  .sensor  coatu.g  layer  inter¬ 
faces;  this  leads  to  a  potential  solution  with  a  completely  different  form.  Since  the  sensor  coating 
modelled  herein  (the  CHEMFET’s  phthalocyanine  polymer)  has  a  fairly  low  transition  frequency 
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(below  1000  Hz),  the  concerns  at  the  lower  frequencies  should  not  represent  a  significant  limitation 
for  the  model. 

5.2  Significance  of  the  Thesis 

While  particular  values  for  the  interdigitated  electrode  structure  have  been  calculated  and 
graphically  displayed,  the  main  benefit  of  this  research  is  in  the  ilexibility  of  the  model  to  handle  a 
variety  of  model  parameters  within  the  investigated  limits  of  the  analysis.  In  this  light,  the  model 
developed  in  this  thesis  becomes  a  powerful  tool  to  further  investigate  interdigitated  electrode  struc¬ 
tures  and  their  use  in  integrated  circuit  sensors.  As  both  Lee  [14]  and  Wiseman  [18]  recommend, 
further  research  in  the  integrated  circuit  sensors  requires  a  rigorous  analysis  of  the  interdigitated 
electrode  structure.  The  model  developed  in  this  thesis  is  intended  to  help  satisfy  this  need. 

8.3  Recommendations  for  Future  Research 

The  ultimate  objective  of  the  integrated  sensor  is  to  detect  electrical  variations  caused  by  an 
external  event.  Therefore,  any  variations  due  to  internal  events  should  be  minimized,  or  negated  if 
possible.  The  depletion  layer  capacitance  (discussed  in  Chapter  III)  caused  by  a  metal-insulatcr- 
semiconductor  contact  presents  an  undesirable,  frequency-dependent  internal  variation  which  coulc' 
hinder  the  sensor’s  fundamental  objective.  Therefore,  this  unwanted  capacitance  should  be  elim¬ 
inated  by  introducing  a  suitable  ground  plane  above  the  semiconductor  substrate.  A  layer  of 
polysilicon,  or  other  appropriate  material,  can  establish  such  a  ground  plane  to  match  the  simpli¬ 
fied  model  developed  in  this  thesis.  Hence,  an  integrated  circuit  sensor  using  the  interdigitated 
electrode  sensor  structure  should  incorporate  this  additional  layer  below  the  structure’s  insulating 
layer. 

The  results  presented  in  this  thesis  also  illustrate  the  importance  of  the  sensor  coating’s 
permittivity  and  conductivity  on  the  overall  results.  Unfortunately,  these  properties  for  the  sensor 
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coating  used  in  the  CHEMFET  are  not  well  established.  While  the  electrical  properties  are  variable, 
the  range  of  this  variation  should  be  determined  to  some  resonable  degree  of  accuracy.  Only  then, 
will  an  acceptable  model  of  the  device  be  available. 

In  addition,  this  research  considered  a  particular  configuration  of  the  interdigitated  electrode 
structure.  Specifically,  the  thin  sensor  coating  layer  was  limited  to  the  region  above  the  electrodes. 
To  further  model  the  CHEMFET,  a  multi-layered  combination  of  P-glass,  polymer,  and  air  should 
be  considered  between  the  electrodes.  In  fact,  different  configurations  of  all  the  layers  within  the 
structure  should  be  explored.  The  solution  technique  employed  in  this  thesis  should  be  readily 
adaptable  to  any  rectangular  geometry.  Alternate  methods  may  have  to  be  used  for  more  complex 
geometries,  such  .as  the  finite  difference  method  used  by  Lee  [14]. 

Finally,  the  unexpected  behavior  of  the  model  element  values  at  low  frequencies,  when  com¬ 
plex  permittivity  is  considered,  should  be  further  investigated.  This  thesis  attempted  to  extend 
the  high  frequency  solution  to  much  lower  frequencies.  Another  possible  approach  would  be  to 
determine  the  potential  everywhere  in  the  periodic  cell  at  DC,  and  then  extend  this  solution, 
up  in  frequency,  to  the  sensor  coating’s  transition  frequency.  For  the  DC  solution,  the  potential 
distribution  is  determined  first  within  the  sensor  coating  layer,  as  discussed  in  Chapter  V.  The 
potentials  in  the  other  layers  are  determined  using  the  calculated  boundary  potential  at  the  sensor 
coating  interfaces  as  boundary  conditions  for  determining  the  potential  in  the  dielectric  layers.  At 
frequencies  below  the  sensor  coating’s  transition  frequency,  the  potential  within  the  cell  can  be 
considered  as  a  perturbation  of  the  DC  solution,  which  will  require  a  new  analysis  and  solution  for 
the  perturbation  term.  This  approach  should  easily  model  the  behavior  at  the  lower  frequencies, 
and,  together  with  the  results  of  this  thesis,  provide  a  complete  frequency-dependent  model  of  the 
interdigitated  electrode  sensor  structure. 
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Appendix  A.  Detailed  Capacitive  Potential  Solution 


To  determine  the  capacitive  elements  of  the  model,  the  potential  distribution  within  the 
simplified  periodic  cell  is  determined  for  the  case  when  the  sensor  coating  is  considered  a  dielectric 
(high  frequency  limit  on  the  boundary  conditions  for  the  sensor  coating).  First,  a  Fourier  series 
representation  of  the  potential  functions  of  each  layer  arc  found  using  unknown  constants.  These 
potential  functions  satisfy  all  boundary  conditions  except  those  at  the  interfaces  between  layers. 
Then,  the  Fourier  coefficients  are  determined  by  enforcing  the  boundary  conditions  between  the 
layers  within  the  cell.  A  diagram  of  the  simplified  periodic  cell  is  depicted  in  Figure  A.l. 


■ 


Figure  A.l.  Problem  diagram  for  the  detailed  capacitive  solution 


A-1 


A.  I  Potential  in  Layer  2 


In  layer  2,  the  voltage  potential,  <^2i  satisfies  (see  Figure  A.l) 


5</>2  A  s  +  1 
2 


4>2{x,-d)  = 


For  <j)2,  we  seek  a  solution  of  the  form 


4>2  =  'i2ix:)^2iy] 


Using  Equation  A. 4  ii  Equation  A.l  yields 


'F  4*  +  =  0 


where  'i'"  denotes  and  denotes  Dividing  by  '!'$  (to  separate  the  variables) 

yields 


(S/» 

Xjf  ~  $  “  ~  2 


where  A2  is  the  separation  constant  or  eigenvalue  in  layer  2.  First,  let  us  resolve  ^(z) 


<?"  +  =  0 


To  simplify  the  formulation,  let 


X  =  X  + 


then  'I'(z)  =  ^(z),  and  d'^'^/dx^  =  d^^/dx^.  Thus, 


+  A?$  =  0 
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or,  we  may  write 


4*  =  cos(A2x)  +  B  sin(  A2X) 

Taking  the  first  derivative  yields 

^ '  =  — ylA2  sin(A2x)  +  BX2  cos(A2x) 

Enforcing  Equation  A. 2  on  Equation  A.  11  yields 

^'(0)  =  flA2  =  0 

or  B  =  0.  Also  from  Equation  A. 2 

$'(s  +  1)  =  -AA2  sin  [A2(s  +  1)]  =  0 

This  condition  forces  the  eigenvalue  A2  to  be 

A2  =  ,  for  n  =  0, 1, 2, ... 

s  +  1 

Thus, 

't'nCx)  =  A„  cos  ( 

or,  using  the  original  x  variable,  the  eigenfunction  is 
^’„(x)  =  AnCos  ^7:^  +  0 

For  the  equations  can  be  simplified  by  letting 

y  =  y  +  d 

Consequently,  <I>(j/)  =  $(j/),  and  =  d^^/dy^.  Then  from  Equation  A. 6 


(A.  10) 


(A.ll) 


(A.12) 


(A.13) 


(A.M) 


(A.15) 


(A.16) 


(A.17) 


<I»o  =  0  for  n  =  0 

=  0  for  n  ^  0 


(A.18) 

(A.19) 
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or 


$0  = 

Coy  +  Do 

(A,20) 

= 

C„  cosh(A2y)  +  D„  sinh{A2y) 

(A. 21) 

Enforcing  Equation  A. 3  yields  (y  =  —d  implies  that  y  =  0) 

io(0) 

=  Do 

(A.22) 

$n{0) 

=  c„ 

(A.23) 

or,  Do  =  Cn 

=  0.  Thus, 

$o(y0 

=  Coy 

(A.24) 

$n(yO 

=  D„sinh(Ay) 

(A.25) 

or 

' 

Co(y  +  d)  for  n  =  0 

'**(y)  =  (A. 26) 

D„  sinh[A2(t/  +  d)]  for  n  ^  0 

Using  Equations  A.  16  and  A. 26  in  A. 4  yields  the  final  form  of  the  potential  in  layer  2  as 

OO 

<i>2  =  ao(y  +  d)  +  ^  o„  sinh[A2(y  +  d)]d’^(a:)  (A. 27) 

n  =  l 

where 


ao 

—  Aq  ■  Co 

(A. 28) 

=  An  -  Dn 

(A.29) 

A2 

nir 

s  +  1 

(A.30) 

(  ^  lA' 

=  cos  nir  - -  +  - 

\  s  +  1  2  / 

(A.31) 
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The  superscript  in  (f>n{y)  denotes  layer  2  instead  of  the  conventional  squared  notation. 


A. 2  Potential  in  Layer  3 

Layer  3  is  located  between  the  two  fingers  with  (0  <  t/  <  /i)  and  1— s/2  <  x  <  s/2).  Referring 
to  Figure  A.l,  the  potential  in  layer  3,  03,  satisfies 

V^03  =  0  (A. 32) 

<^3(-|,y)  =  V^i  (A. 33) 

<^3(+|,y)  =  V7  (A.34) 

For  the  potential,  we  seek  a  solution  of  the  form 


y)  =  V3{x)  +  W3(x,  y) 


(A. 35) 


where  t)3(x)  is  selected  to  satisfy  Laplace’s  equation  and  the  left  and  right  boundary  conditions. 
Let  us  consider  the  following  function  for  t;3(x) 


1^3(1;) 


V^2-V^i  ,  ^.  +  Ri 

- ^ — 


(A. 36) 


This  function  satisfies  Equations  A. 32,  A. 33  and  A.34.  Thus,  W3(x,y)  must  satisfy 


V^W3 

=  0 

(A.37) 

w3(-s/2,y) 

=  0 

(A. 38) 

W3(+s/2,y) 

=  0 

(A. 39) 

To  solve  W3(x,  y),  we  seek  a  solution  of  the  form 

W3(x,y)  = '*'3(x)^3(y)  (A. 40) 

Using  Equation  A. 40  in  Equation  A. 37  yields 
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(A.41) 


where  and  <f"  =  Thus, 

’*’"  _  _  \2 
$  “  $  ~  3 


(A.12) 


where  A3  is  the  separation  constant  or  eigenvalue  for  layer  3.  The  analysis  can  be  simplified  by 
letting 


Then  ^(x)  =  'I'(x),  and  /dx^  =  5^'f'/5x^,  so 


(A.43) 


$"  +  =  0 


(A.44) 


or 

$(x)  =  A  cos(A3x)  +  fl  sin(A3x)  (A. 45) 

From  Equation  A. 38 

$(0)  =  A  =  0  (AA6) 

From  Equation  A. 39 

^(s)  =  fl  sin(A3s)  =  0  (A. 47) 


or 


A3  =  — ,  for  n  =  0, 1,2, ... 
s 

Hence  the  eigenfunction  for  layer  3  is 


^3,(?) 

'f3,(x) 


Bn  Sin( - ) 

s 


B„  sin 


n?r 


{A.48) 


(A. 49) 
(A. 50) 
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F  from  Equation  A. 42 


<1."  -  =  0 

Since  we  expect  some  symmetry  about  y  =  let 

h 

y  =  y-2 

Then  <^(y)  =  ^(y),  and  =  d^^/dy^.  Thus, 

^3(y)  =  Cn  cosh(A3y)  +  Dn  sinh(A3y) 

Equivalently  (since  ^'33  =  0,  was  not  sought), 

‘^3(y)  =  Cn  cosh  ■^3^y-^^  +£>„sinh  -^3  ^y  -  ^  j 
Using  Equations  A. 50  and  A. 54  in  A. 40  yields 

•^3(-c,y)  =  ^  cosh  +c„sinh  -^3  ^y  -  | 

where 

tPl(x)  =  sin  nn  i  j 

6fi  ~  AnCfi 

Cfi  —  An  Dn 


(A.ol) 

(A, 52) 

(A. 53) 

(A.54) 

(A.55) 

(A. 56) 
(A.57) 
(A. 58) 


The  superscript  in  denotes  layer  3.  Finally,  from  Equation  A. 35,  the  potential  in  layer  3  is 

<?^3('E,y)  =  «^3(i)  +  W3(x,y)  (A. 59) 

with  t»3(x)  and  W3(z,y)  defined  in  Equations  A. 36  and  A.55,  respectively. 


A-7 


A. 3  Potentials  m  Layers  4  and  5 


Layer  4  is  located  immediately  above  the  electrodes  (h  <  y  <  /),  while  layer  5  is  above  layer 
4  (/  <  j/  <  oo).  In  each  of  these  layers,  the  potential,  (f>4  or  <^5,  satisfies  (see  Figure  A.l) 


=  0  forJfc  =  4,5 

^(±1^,;,)  =  0  rotl-  =  4,5 


(A. 60) 
(A,61) 


The  potential  also  satisfies  the  boundary  conditions  developed  in  Chapter  III.  From  Equa¬ 
tions  3.57,  3.66,  and  3.72 


lim  |05(z,y)|  < 

f  — OO 

OO 

(A.62) 

<i>A{x,l)  = 

4>s{x,l) 

(A.63) 

II 

d<i>s  1  ,, 

(A.64) 

Potential  in  Layer  5.  For  ^5,  we  seek  a  solution  of  the  form 
<t>s  =  'lf5{x)^5iy) 

Using  Equation  A. 65  in  Equation  A. 60  yields 

=  0 

where  denotes  and  <I>"  denotes  .  Thus, 

^  ^  _  i2 

^  ~  $ 

First  'I’(x),  from  the  previous  equation 

=  0 


(A.65) 


(A. 66) 


(A.67) 


(A. 68) 
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where  A5  is  the  separation  constant  or  eigenvalue  for  layer  5.  Since  this  situation  is  identical  to 
in  layer  2,  including  the  boundary  conditions,  the  results  of  'Pt,  can  be  used.  From  Equation  A.  16 


=  An  cos 


X  1 

nir  I  - +  - 

s+  1  2 


and 


nir 


As  =  - ,  for  n  =  0,  1 , 2 _ 

s  +  1 


Now  for  from  Equation  A.C 


<!►"  -  A^<I>  =  0 


=  0  for  n  =  0 


0"  _  A^$„  =  0  for  n  ^  0 


Thus, 


‘^*o(j/)  =  C'oy  +  Do 

^Ay)  =  Cn  exp(A5y)  +  Dn  expf-Asy) 

Enforcing  Equation  A. 62  yields 


lim  |Coj/  +  Dol  < 

y  — *00 

lim  iC„exp(A5yi  +  E:nexp(-A5y)|  <  00 

y  — •  CO 


Hence  Cq  =  Cn  —  0,  and 


‘I>n(y)  =  { 


Do 


Dn  exp(-A5y) 


for  n  =  0 
for  n  ^  0 


(A. 69) 


(A. 70) 


(A, 71) 


(A, 72) 


(A.73) 


(A.  7-1) 


(A. 7.5) 


(A. 76) 


(A. 77) 


(A.78) 
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From  Equation  A. 65 


J/)  = (A. 79) 

or 

oo 

(isl-i’.!/)  =  Do  +  ^  'I'n(i:)D„  exp(-A5y)  (A, 80) 

fl  =  1 

since  'I'5„(i^)  =  -lo  and  the  constants  were  combined  into  Dq. 

Potential  in  Layer  1,.  For  d>4(r,i/),  we  again  seek  a  solution  of  the  form 
=  'l'4(-r)<I>4(i/)  (A. 81) 

('sing  Equation  A. 81  in  Equation  A. 60  yields 

=  0  (A. 82) 


or 


vjt"  4)" 

v|r  “  "  <J)  ~ 


(A. 83) 


where  A4  is  the  separation  constant  or  eigenvalue  for  layer  4.  First  determine  'F4(r),  from  the  above 
equation 


♦"4-  A^'t'  =  0 


(A. 84) 


Since  this  situation  is  identical  to  'Fsfr),  including  the  boundary  conditions,  then 
'I'4,  =  'I's,  =  A„  cos  j^riTT  ^  0 

and 


(A  85) 


A4  =  A5  =  - ,  for  n  =  0,  1.2,... 

s  +  1 


(A. 86) 
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Now  for  <I>4,  let 


y-y-i  (A. 87) 

which  leads  to  <^(y)  =  ^(y),  and  d~^/d^  =  d^^/dy^.  From  Equation  A. 83 

$"-A-$  =  0  (A. 88) 


'I>Q  =  0  for  n  =  0 


C  - 


0  for  n  0 


(A. 89) 
(A, 90) 


Hence 


‘^n(y)  =  { 


F'oy  +  Eo  for  n  =  0 

En  cosh(Ay)  +  Fn  sinh(A(/)  for  n  ^  0 


(A.91) 


and 


y)  -  Foy  +  Eq  +  ^  'Fn  [En  cosh(AijO  +  E„  sinh(Ay)] 


A.  92) 


n=0 


Applying  Boundary  Data.  Using  Equations  A. 92  and  A. 80  in  A. 63  yields 


Eo  +  ^  'Fn(i)E„  exp(-A/)  =  Eo  +  ^  'F„(x)En 


(A  93) 


n  =  l 


n=0 


Thus, 


Eo  =  Do  (A. 9-1) 

En  =  D„exp(-M)  (A  9.3) 


Using  Equations  A. 92  and  A. 80  in  A. 64  yields 
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fs  [-AD„]exp(-A/)  =  (4 


n  =  1 


Fo  +  ^  'I'n(-r)AF„ 


n=:0 


Thus, 


Fo  =  0 


Fn  =  --D„exp(-A/) 

f4 


Hence, 


?4(-r,  ^  -  Do  +  ^'f/nD„  exp(-A/) 

n  =  l 

or  the  potential  in  layer  4  can  be  expressed  as 

00 

<>4(x,  y)  =  do  +  Y^  d„Tn{y)^ni^) 

where 


cosh(At/) - ^sinh(A^ 


n  =  l 


K{y) 

dn 

do 


=  exp(-A4/)  \  cosh  [A4(/  -  y)]  + 


—  AnDn 


f5 


sinh[A4(/-y)]| 


Do 


cos 


X  I 

nrr  1 - f-  - 

s  +  1  2 


(A. 96) 


(A. 97) 


(A, 98) 


(A. 99) 


(A. 100) 


(A. 101) 
(A. 102) 
(A. 103) 
(A. 104) 


Finally,  the  potential  in  layer  5  is 


00 

<i>^{x,  y)  =  do  +  Y^dn  exp(-A5y)^^(x) 

n  =  l 


where 


f  ^ 

COS 

nn  - 

+  - 

L  v«  +  i 

2yJ 

(A. 105) 


(A.  106) 


The  superscripts  on  ip^  and  denote  the  material  layer  of  interest  instead  of  the  conventional 
power  operator  notation. 
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.4.^  Boundary  Between  Layers  2  and  3 


Summary  of  Potential  Functions.  From  Equations  A. 27  A. 36,  A. 55  and  A. 54 


4>2(x,y)  =  ao(y  +  d)  +  ^a„sinh[A2(y  +  d)]i/'2(x) 

n  =  l 
oo 

Oo  +  ^  cosh[A2(y  +  d)il<l(x) 

n  =  1 

l-'3(^)  +  jftn  cosh  +CnSinh  |  t  nl 


502  .  . 


03(x,y)  = 


(A.  107) 
(A. 108) 
x)  (A. 109) 


502 

-5— (-c.l/ 
5y 


=  ^A3|6nsinh  -'a  ^y  -  +  <•„  cosh  . 


>a  )  y  -  - 


V'^(r)  (A. 110) 


where 


V3{x) 

A2 

-^3 


COS 


X  1 

nir  I  - -  +  - 

s  -|-  1  2 


,  X  I 


sm 

'^2-  V^l 
s 
nir 


X  + 


V2  + 


s  +  1 

nir 

s 


(A. Ill) 
(A. 112) 
(A. 113) 
(A. 114) 
(A. 115) 


Boundary  Conditions.  The  boundary  conditions  for  the  interfaces  between  layers 
were  developed  in  Chapter  III.  From  Equations  3.58,  3.60,  3.64,  and  3.70 


02(1,0)  =  f{x) 


(A. 116) 


502.  503 


(A. 117) 


where 
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/(x)  = 


V^I,  for  -i±l  <  X  < -f 

03(x,  0),  for  -5  <  X  <  I 

^2, 


(A. 118) 


for  I  <  a:  < 


Using  Equations  A. 107  and  A. 109  in  A. 116  gives 


aod+  ^ 


(A. 119) 


n  =  l 


k'here 


0'  =  sinh 


mrd 
s  +  1 


(A. 120) 


To  determine  the  a„  terms,  Equation  A.  119  is  “sifted”  by  multiplying  each  side  by  and  inte¬ 
grating  with  respect  to  x  over  the  interval  Accordingly, 


aod(l,l)  =  (/(x),l) 

00 

=  {/(x),V’m) 


(A. 121) 
(A. 122) 


nstl 


where  “1”  denotes  the  function,  tpQ,  which  has  constant  unit  value  over  the  interval  and 

tpQ  =  1  for  m  =  0 

{g(x),h(x))  =  J  ^g{x)  h{x)dx 


(A. 123) 
(A  124) 


(/'“  are  cosine  eigenfunctions  which  ace  orthogonal  to  the  constant  function.  Therefore,  only  the 
rn  =  0  term  is  used  for  sifting  the  constant  term.  First,  evaluating  (1,1)  yields 


/lil 

<fx  =  (s  -t-  1) 


Evaluating  (/(x),  1)  yields 


(A. 125) 


(/(■f),  1)  =  (Vi ,  1)-  +  {<^3(x,  0),  1)0  +  {U2.  1)+ 


(A. 126) 
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where 


(V'l.l)-  =  Vi 


and 


{V2,1)+  =  V2 1 


(A. 127) 


(A. 128) 


Next, 


{4>3{x,  0),  l)o  =  l)o  +  ^  (V>^,  l)o 


n  =  l 


where 


<*>'  = 


n  = 


cosh 

.  ,  /  nirh 
sinh  — — 

V  2s 

■*  'V2-V, 


dx 


Using  ^  —  7  +  I  yields 


-i: 


sin(nirr)s  dr 


=  -::^[<^os(n7rr)]‘ 


nn 


-— [cos(n7r)  -  Ij 
mr 


(A. 129) 


(A. 130) 
(A. 131) 


(A. 132) 
(A. 133) 


(A. 134) 


Using  Equations  A. 125,  A, 126,  A. 127,  A. 128,  A. 129,  A, 132,  A. 134  in  Equation  A. 121  yields 

aorf(s+l)  =  ^  +  ^  + £(!/,  + 1/2)+^  ^  [!-(-!)"]  (A.I35) 

n  =  l 
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or 


aod(s  +  1)  -  ^  [6„$=  -  c„n]  ^  [1  -  (-1)"]  =  ^  ( V,  +  V,) 
;  nir  I 

n  =  l 

Next,  evaluating  yields 


To  simply  the  algebra,  let 

X  1 
— T  +  :^ 


then 


X  =  (s  +  l)r  —  t 
dx  =  (s  +  1)  dr 


Using  Equations  A.  138  and  A.  140  in  A.  137  yields 


=  (s+l)  cos(n5rr)cos(rrnrr)(ir 

Jo 

s+l 

—  -  /  cos(na) cos(ma) da 

Jo 

where  a  =  wt  and  a  =  Trdr.  For  integers  m  and  n 
/  cos(nx)  ■  cos(mr)  dx  = 

Jo 


0,  for  n  ^  m 
for  n  =  m 


Thus, 


0,  for  n  ^  m 


for  n  = 


m 


Evaluating  {f{x),4>m)  yields 


(A. 136) 


(A. 137) 


(A. 138) 


(A. 139) 
(A. 140) 


(A. 141) 
(A. 142) 


(A. 143) 


(A. 144) 
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=  {ViAV)-  +('^3.V?„)'’  +  (V'2,^m>+  (A.  145) 

where 


=  V,  f  ^  rkldx 

(A. 146) 

(<^3,V'm>0  =  J  ^4)3{x,0)xl)l,dx 

(A. 147) 

3 

{V,,rPl)+  =  i>ldx 

(A. 148) 

For  (Vi.V'm)-.  using  Equations  A. Ill,  A. 138  and  A. 140  in  A. 146  yields 
('^i.  V'm)-  =  (*  +  l)Vi  /  cos(m7rr)dr 
where  r'  =  -  +  i  =  0  and  r/  =  -  +  i  =  5^.  Integration  yields 


=Ki(s  +  l) 


- sin(m7rr) 

mTT 


Jo 


or 


mTT 


For  (V2>V'm)  +  .  using  Equations  A. Ill,  A. 138  and  A. 140  in  A. 148  yields 
V’m)+ =  («  +  1)V2  j  ^  cos(m7rr)(ir 

where  2  ~  i»+2  ~  ^2^J4T  +  j  =  1-  Integration  yields 

1 


(v^2.<r  =  v^2(s+i) 


mir 


sin(mT7') 


2  j  +  l  ^  2 
1 


or 


V2{s+l)  . 

- sin 

mir 


mir 


f  — ) 

\2s+2j 


Since 


(A. 149) 


(A. 150) 


(A. 151) 


(A. 152) 


(A. 153) 


(A. 154) 
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(A. 155) 

(A. 156) 

(A. 157) 

(A. 158) 

(A. 159) 

(A. 160) 
(A. 161) 

(A. 162) 
( A.  i63) 

(A. 164) 
(A. 165) 
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Using  Equations  A. 144,  A. 145,  A. 151,  A. 158,  A. 159,  A. 163  in  A, 122  yields 


dfri  ^ 


c 

m 


S  +  1 
2 


where 


=  +  -  CuK]  fL 

n  =  l 


(A.  166) 


(A. 167) 


or 


<2rn  ^ 


C 

m 


S  +  1 
2 


n=  1 


(A. 168) 


For  the  jump  condition,  using  Equation  A. 108,  A. 110,  A. 114,  and  A. 115  in  A. 117  yields 


f2 


ao  +  ^ 


n  =  l 


riTT 
5  4-1 


+  C„«I>' )  tp] 


n  =  l 


where 


(A. 169) 


(A,  170) 


Here,  sifting  occurs  with  ip^,  because  the  interval  over  which  the  jump  condition  applies  is 
-s/2  <  X  <  s/2.  Thus, 


f2ao(l,  'Pm  )  +  E  =  C3  ^  ^  )  (0^0^) 


fl  =  l 


n  =  l 


(A.171) 


which  leaves  (1,0^),  (0^,0^)  and  (0^,0^)  to  be  solved.  First,  (1,0^)  can  be  written  as 


(A.  172) 


Using  r  =  j  +  i  and  sdr  =  dx  yields 


(1,0^)  =  si  sin(mTr)dr  = 

Jo 


cos{mirr) 


mir 


=  — [- cos(m?r)  +  1]  =  —  [1  -  (-1)'' 

mJT  rrjT 


(A,.  73) 


Evaluating  yields  (see  Equation  A. 163) 


\  (  ^  iM  •  r  iM 

=  /  cos  TIT  - H —  1  ■  sin  mT  —  +  - 

J-sn  V*  +  i  2)\  [  \s  2)\ 


(A. 174) 
(A. 175) 


Next 


h  {~s + 0]  H  ^  0] 


Using  r  =  j  +  i  and  sdr  =  dz  yields 


(A. 176) 


(V'niV’m)  ^  /  sin(nTr)  •  sin(mTr) dr 

Jo 

s 

=  —  /  sin(na)  •  sin(ma)  dot 

do 


(A. 177) 


where  a  =  tt.  For  integers  m  and  n 


sin(nz)  •  sin(mx)  dx  = 


0,  for  n  ^4  m 
f ,  for  n  =  m 


Thus, 


0,  for  n  ^  m 
f ,  for  n  =  m 


Using  Equations  A. 173,  A. 175,  and  A. 179  in  A. 169  yields 


(A. 178) 


(A. 179) 


[1  -  (-1)-]  ao  +  C2  ^ 


(A. 180) 


Summary  of  Equations  for  Layer  2-3  Interface.  The  results  of  matching  the 
boundary  conditions  along  the  interface  between  layers  2  and  3  provided  the  following  equations 
that  must  be  solved  to  determine  the  potential  coefficients 


i^(l/i  +  l/2)  =  aod(s  +  l)- V— (6„<I>'  -c„'Ir')[l-(-l)"] 

A  JIT 


(A. 181) 


(A.  182) 


.4.5  Boundary  Between  Layers  3  and  4 

Summary  of  Potential  Functions.  From  Equations  A. 100  and  A. 109 


00 


=  do  +  J2'^nK{y)Ki^) 

fi  —  1 

(A. 184) 

00 

=  -Y.X,dMy)ri^tix) 

n  =  l 

(A. 185) 

=  1^3(1)  +  ^  |6r.  cosh  +c„sinh  A3  |  r/-: 

a(x)  (A.186) 

=  ^AajAnsinh  A3  +c„cosh  ^3  ^y  - 

(A. 187) 

where 


Kiy) 

Kiy) 

A4 

-^3 


(A, 188) 
(A. 189) 
(A, 190) 
(A. 191) 
(A. 192) 
(A. 193) 
(A. 194) 
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Boundary  Conditions.  The  boundary  conditions  for  the  interfaces  between  the 
I-'.ye-s  were  developed  in  Chapter  III.  From  Equations  3.59,  3. Cl,  3.65,  ^iiid  3.71 

(j)4{x,h)  =  g{x)  (A. 195) 


d<p4 

>n)  =  f3^{x,n) 
oy 

(A. 196) 

where 

' 

V'l,  for  —  —  5 

</(j;)  =  < 

<p3{x,  h),  for  -  |  <  r  <  1 

(A. 197) 

V2,  for  5  <  a;  <  ^ 

Using  Equations  A. 184  and  A. 186  in  A. 195  gives 

OO 

do  +  Yl  =  g{x) 

(A. 198) 

n  =  l 


As  at  the  boundary  between  layers  2  and  3,  the  d„  terms  are  determined  by  “sifting”  Equa¬ 
tion  A.  198.  Each  side  of  the  equation  is  multiplied  by  and  then  integrated  with  respect  to  x 
over  the  interval  — ^  x  <  Since  =  ip^,  the  results  of  Equations  A. 125,  A. 126,  A. 127, 
A. 128,  A. 129,  A. 132,  and  A. 134  can  be  used  in  Equation  A. 198  to  obtain 

‘^0(5-1- 1)  =  q-Kz)  4- f;[6„<I»'+c„n]  —  [l -(-!)"]  (A.iuu; 

z  ^  ^  riTT 

n=l 

=  Im  +  f2i^nK+CnK]lln  (A.200) 

n  =  l 

For  the  jump  condition.  Equations  A.  185  and  A.  187  can  be  substituted  into  Equation  A.  19b  to 
obtain 

oo  oo 

-  (A. 201) 

,  S  T*  1  ,5 

n  =  l  fi=l 
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Here,  as  in  the  lower  boundary  between  layers  2  and  3,  sifting  is  accomplished  using  the  inner 
eigenfunction  over  the  interval  —s/2  <  jr  <  s/2.  Since  =  i'n,  The  results  of  Equations  .4.175 
and  A. 179  can  be  used  to  obtain 


riTT 

5  4-  1 


dnKWiL 


mir  s  , 

s  I 


(A, 202) 


Summary  of  Equations  for  Layer  5-/  Interface.  The  results  of  matching  bound¬ 
ary  conditions  at  the  interface  between  layers  3  and  4  provides  the  following  equations,  which  the 
potential  coefficients  must  satisfy: 


s  -t-  1 


(V,  +  V,)  =  dois  +  l)-f2~if>nK+CnK)[^-(-^n 


n  =  l 


Im  = 


,  ,  OO 

+  CnK)I, 


3 

mn 


n  =  l 


0  = 


fi=i 


A.  6  Evaluation  of  I}, 


From  Equation  A.  164 


(A.203) 
(A. 204) 
(A, 205) 


(A. 206) 


Using  T  =  -f-  i,  and  da:  =  (s  -I-  l)dr  yields 


j  (*  +  1)^  cos(m»T)dr 


(V'2-Ui)(s+1)2 


[/: 


7-cos(mTr)  dr 


-u: 


cos(mnT)  dr 


(A. 207) 


Since 


/ 


x(cos  ax)  dx  =  —  cos  ax  -i —  sm  ax 


(A. 208) 
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1/2  +  Vi  fi 


/  COS  - -I —  1 

V^+i  2yJ 


To  simplify  the  algebra,  let  r  =  +  i.  Then  [s  +  1)  dr  =  dx  and 


r2  ''2  +  n 

Im  =  - ^ - J  cos{mwT){s  +  l)dr 

where  n  =  2j+2'  ~  §7+1  Integration  yields 

•>  1  r  1  V‘ 

=  o(^2  +  Vi)(s  +  1)  - sin(mirr) 

z  mTT 

*“  Jr, 

(V2  +  Vi)(s+1)  f  .  /  2s+l\ 

=  - - -  sin  mir -  —  i 

2m7r  \  2s +  2/ 


(V2  +  Vi)(s+1)  [  .  /  2s +1 

- - -  sin  mir - 

2m7r  \  2s +  2 


2s +  2 


,2s  +  1 
2s  +  2 


=  -(-ir  sin 


2s  +  2 . 


-•1..?  Evaluation  of 


From  Equation  A.  163 


/L  =  /_%in  [n.  +  i)]  cos  [m.  +  i)]  dx 


To  simplify  the  algebra,  let 


X  1 

1-  =  -  +  - 
s  2 


then  X  =  ST  —  and  dx  =  sdr.  Hence, 


=  /  sin(n5rr)cos  mir  [ - - —  +  -  |  sdr 

Jo  V-s+l  2s +  2  2/_ 

f'  ■  ,  .  /2sr+n]  , 

=  s  /  sin(n7rr)cos  mir  I  -  dr 

Jo  V  2s  +  2  y 


(A. 211) 


(A  .215) 


(A.216) 


(A.217) 


(A. 218) 


(A. 219) 


(A.  220) 


(A. 221) 
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To  further  simplify  the  formulation,  let  fh  =  If  n  =  fh,  then 
1  f’' 

Imn  =  -  cos(na)sin(na)cfor  (A. 226) 


Since 


cos  ax  ■  sin  ax  dx  =  0 


(A.227) 


then 


I 


3a 

mn 


=  0,  for  n  = 


ms 
5  +  1 


(A, 228) 


For  n  in 


J 


(cos  az)  (sin /?x)  dx  =  -- 


cos(/i  —  a)x 
2(/?-a) 


cos(/?  +  a)x 

2[P  +  n) 


then 


(A. 229) 


I 


3a 

mn 


cos  (ma)  sin  (na)  da 
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1  cos(n  — m)a  cos(n  +  m)ol" 

TT  2(n  —  m)  2(n  +  m)  g 

_  If  cos[(n  — m)7r]  cos[(n  +  m)T]  1  1 

TT  \  2(n  —  fh)  2(n  +  m)  2  (n  —  rn)  2  (n  +  m) 


Since  cos(a  +  /?)  =  cos  a  cos  /?  —  sin  a  sin  0,  then 


n  +  ’7i  +  n  —  m  1 


—  [cos(n7r)  cos(mT)  +  sin(rnr)  sin  (rnTr)] 


2(n2_m2)  . '  '  ' 

- - ^  [cos(n7r)  cos(mir)  +  sin(n)r)  sin(m7r)] 

2(n  +  m) 


Using  sin(n7r)  =  0  and  cos(n7r)  =  (—1)"  yields 


C,  =  -j — fi  -  (-1)” (^)l .  f»'  ^  #  a 

For  from  Equation  A. 224 


=  /„ 


.  mirsT\ 

sin  -  1  sin(n7rr)ar 

\s+l  J 


Again,  for  simplification,  let  a  =  nr.  Then  dr  =  —,  and 


Using  m  =  yields 


1  /•' 

=  —  /  sin  (ma)sin(na)  do 

Jo 


If  n  =  fh.  then 


=  —  f  s\n^ (na)  da  =  ^  f  [l  —  cos(2na)]  da 

TT  Jo  2ir  Jq 

=  —  [q - sin(2nQ)]  =  —  (;r  —  0  —  0  +  0) 


Thus, 
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(A. 245) 


(A. 246) 


(A.247) 


(A. 248) 


A. 9  Summary  of  Linear  Equations 

To  determine  the  potential  coefficients,  the  results  of  matching  the  boundary  conditions  at 
the  interfaces  between  layers  2-3  and  3-4  can  be  consolidated  by  combining  the  resulting  equations 
at  each  interface.  This  will  establish  a  final  set  of  equations  that  the  potential  coefficients  must 
satisfy  to  fulfill  all  boundary  conditions.  Adding  Equations  A. 181  and  A. 203  yields 
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(A. 249) 


aod(s  +  1)  -  V  — [1  -  (-1)"]  +  dois  +  1)  =  (s  +  1)  (l/i  +  \/2 

HTT 

n  =  l 

Subtracting  Equation  A. 203  from  A. 181  yields 


aod{s  +  1)  +  y;  —c„K  [I  -  (-1)")  -  do(s  +  1)  =  0 

HTT 

n  =  l 

Adding  Equations  A.  182  and  A. 204  yields 


«me^  =  21. 

n  =  1 

Subtracting  Equation  A. 204  from  A.  182  yields 


.s  +  1 


n  =  l 


^4/A^£±I  _ 


Adding  Equations  A. 183  and  A. 205  yields 


n  =  l 


^  [1  -  (-l)-]  ao  +  J2  +  Yi 

ns  1  —  ’ 

Subtracting  Equation  A. 205  from  A.  183  yields 


^  [1  -  (-I)”|  a,  +  x:  -  Y,  = 


e2n)r 


■w4/lw3  _ 


n  =  l 


A. 10  Glossary  of  Variable  Quantities 


The  following  Quantities  were  used  to  simplify  the  final  expressions: 


©^ 

r' 

Kih) 


sinh 

cosh 


f  mird 

VT+tJ 

mird  \ 

7+1 J 


=  exp  {-Xml)  |cosh  [Am(/  -  ^)]  +  ^  sinh  [Am(/  -  h)] 
=  exp  {—X„l)  |sinh  [A„(/  -/»)]+—  cosh  [A„(/  -  /i)]| 


=  cosh 


nirh 

~27 


(A.  250) 

(A. 251) 

(A. 252) 

0  (A. 253) 

0  (A. 254) 


(A. 255) 
(A. 256) 
(A.257) 
(A. 258) 
(A.259) 
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(A. 260) 


=  sinh 


A„  = 


(^) 


{Vurl^l)-  = 


^(s+1) 


.  f  mn  \ 

sin  ) 

\  2s  +  2/ 


(V^2-^)(s+l) 


“"(5r?2)f-'-'»"i 


(A.261) 


(A. 262) 


(A.263) 


(A.264) 


(A. 265) 


(A.266) 


(A. 267) 


where  the  superscript  c  denotes  the  capacitive  potential  solution.  In  addition,  m  and  n  are  integers. 
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Appendix  B.  Detailed  Charge  Calculations 


To  find  the  total  charge  on  either  microstrip,  the  normal  electric  flux  density  is  integrated 
around  a  contour,  which  encloses  the  microstrip.  A  diagram  of  the  simiplied  periodic  cell  is  depicted 
in  Figure  B.l,  which  also  shows  the  contour  used  for  the  charge  calculations. 


Figure  B.l.  Contour  for  Charge  Integration 


Thus,  the  charge  on  the  left  electrode  is 
Q  —  t  j)  En  dl 

where  7  =  cl  +  c2  +  c3  +  c4,  which  is  shown  in  Figure  B.l.  This  yields 

Q  =  f2  /  {-Ey)dx  +  €2f  Ej:dy  +  (3  f  E^  dy 
J-t^  y_i  Jo 

+€4  /  Ecdy-\-tA  [  Ey{-dx) 

Jh  Jo 


(B.l) 


(B.2) 
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Since  —  —d4)/dx,  and  Ey  =  —d<j>/dy,  then  with  substitution 


(B.3) 


or 


Q  —  Q2  +  Q2  +  Q3  +  Q4  +  Q4 


(B.4) 


where  each  integral  in  order  has  been  assigned  a  corresponding  label.  From  Equations  A. 27,  A. 108, 
A. 109,  A  100,  A. 185,  A. 188,  and  A. 189 


d4>2 

dy 

d<t>2 

dx 

d(j)3 

dx 


d4>A 

dx 


<■"  + 1:  “-nwco.  [»•  +  \) 

n  =  l  '  ' 

n  =  l  '  ' 

Yisl}^]  +  Y  (— 

s  )  v  5 


n  =  l 


{  bn  cosh 

HTT  /  /l\l 

nn  f  h\] 

1 

fx  IN' 

—  y  -  0  +Cn  sinh 

>  cos 

ntr  (  -  +  -  ) 

1 

L «  V  2yj 

L  *  \  2/J 

J 

L  Vs  2/J 

=  -E 


HTT 
5  +  1 


c/„7-„^(y)sin  [n^  (7^7  +  ^) 


(B.5) 

(B,6) 


(B.7) 

(B.8) 

(B.9) 


where 


r;(y) 

=  cosh[A(y  +  d)] 

(B.IO) 

0n(y) 

=  sinh[A(t/ +  d)] 

(B.ll) 

^n(y) 

=  exp(-A/)  (cosh[A(/  -  y))  +  —  sinh(A(/  -  y)j ] 
L  (2  J 

[  (B.12) 

^^(y) 

=  exp(  — A/)  1  sinh[A(/  —  y)]  +  —  cosh[A(/  —  y)]"] 

1  ^2  J 

[  (B.13) 

A 

nir 

s  +  1 

(B.14) 
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First  Ql,  from  Equations  B.3  and  B.4 


Ql 

Ql 

For 


Q2'  —  2^2flo(s  +  1) 


For  Qf 


=  f2  ^OnF^  [sin  (y) -®*"(0)] 

n  =  l  '  ' 


Using  Equations  B.17  and  B.18  in  B.16  yields 
Ql  =  ^f2ao(s  +  l)  +  f2^  a„r'  ^-0sin 

n  =  l  '  ' 

For  Q5 

/•o  00  /  \ 

B-3 


r 

£2  /  aodx 

+^2  E  «nr'  f-0  cos  +  i) 

•'-•  r-  n  =  i  '  '  V  •/  V-’  •  ^  -/ J 


=  Qf+Qf 


(B.15) 

(B.16) 


(B.17) 


(B.18) 

(B.19) 


Evaluating  the  integral  separately  yields 


/  *^niy)<^y 

J-i 


r 
l_ 

s  +  1 


/  sinh[A(y  + 

J-i 


/ 


®n(y)dy  = 


i 


nir 

s  4-  1  [” 
nir  [ 


-4 


cosh[A(y  +  d)] 

n(0)-n  (-^ 


Using  Equation  B.21  in  B.20  yields 


oo  P  /  1 

n=:l  L  \  /  . 


By  letting 


Q2  =  Qt  +  Q^2 


then 


or 


<32  =  ^f2ao(5+ l)  +  €2^^anr'  ^-^sin^^) 
+e2  2a„sin(Y)  r'(O)-r'^-0j 

1  °° 

<32  =  2f2ao(s  +  1)  +  f2  X]  “n  (v) 


n  =  l 


Since 


r' (0)  =  r=  =  cosh 


/  nird  \ 

v^/ 


Now  for  Qs 


(B.20) 


(B.21) 


(B.22) 


(B.23) 


(B.24) 


(B.25) 


(B.26) 
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<33  -  <?3  +  Qz 


(B.27) 


where 


Q3  = 


Q3  = 


-)  dy 


(B.28) 


rh  00 

/  /nr> 

-^3  /  >  - cos  — 

Jo  S  \2  y 

) 

l^nCOsh  '^3^y-^[j 

+  c„sinh  A3  (y-  ^)  1  dy 

(B.29) 

First,  evaluating  yields 


Q3  =  -(V^1-V^2)MS  =  —  (V^i- v^2) 

5  S 


(B.30) 


Now  for  <33 

<33  - 


^  I^COS  ^  |i^<.osh  +c„sinh  -^3  |  (iy 

-^3  E  y  cos  (I^)  ^  |6,sinh  A3  {y  -  ^)]  +  c„  cosh  [A3  (2/  "  ^)]  }[ 
_,3|;cos(^) 

•  (?)  (-?)  (?)  (-?)_ 
-2f3E'^03  ( 

fl  =  l  '  ' 


( nnh\  (  nirh 

+  e„cosh(  — l-c„cosh(-  — 


(B.31) 


Using  Equations  B.30  and  B.31  in  B.27  yields 

<33  =  Y  “  ^^2)  -  2e3  Ecos  (y)  *nsinh 

a  =  l  ^ 


(B.32) 


Now  for  <35 


-  L=A  00  /  \ 

«<  =  ‘'I' ‘‘y 


(B.33) 


Evaluating  the  integral  separately  yields 


LiA  i-Ji 

^niy)dy  =  J  exp(-A/)  |cosh[A(/ -  y)]  +  ^sinh[A(/ -  y)]|  dy 


-exp(-A/)  I  |sinh[A(/  -  y)]  +  —  cosh[A(/  -  y)] 


nw 


/: 


Kiy)dy  = 


riTT 

s  -\-l 

nir 


KW-K 


(h) 


Using  Equation  B.34  in  B.33  yields 

<54  =  f4  ^<i„sin 

n=l  t  \  / . 

Next,  evaluating  QJ  yields 


^4 


cos 


X  1 

nn  1 - f-  - 

5-f-  1  2. 


dx 


i±* 

3 


X  iM 

riTT  1  — — r  +  - 


.J  0/1 


nsl 

oo 


,+  1''""^  2 


I'l  s+i  .  r  ( 

- sin  mr 

J  nir  [  \ 


.5  +  1  2 


X  1 
+  r 


s  +  1  2 


n  =  l 


=  E  ( S- )  h  (t)  ■ 


sin  (y) 


By  letting 


_i±i 


(B.34) 


(B.35) 


(B.36) 


Q4  =  QU  <34 


(B.37: 


then  substituting  Equations  B.35  and  B.3b  in  B.37  yields 


Q4  =  f4^d„sin(Y)  -■^n 

n  =  l  ^ 


+  f4 


n  =  l 


(B.38) 
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Hence, 


OO 

X^rfnsm  (y)  (B.39) 

n  =  l 

Total  Charge.  The  total  charge  enclosed  on  the  left  electrode  (see  Figure  13.1) 
can  be  determined  by  substituting  Equations  B.25,  B.32,  and  B.39  into  Equation  B.4,  which  yields 


1  °°  L 

Q  =  2^2ao(s  +  1)  +  f2^ansin  r'  + -^(V'l  -  Vo) 

n  =  i 

-2f3  £  cos  b„  sinh  f  +  (f4  £  sin 

n  =  l  '  n=l 


(B.40) 
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Appendix  C.  Detailed  Resistive  Potential  Solution 


To  determine  the  resistive  element  of  the  model,  the  potential  distribution  within  the  si  u- 
sor  coating  layer  is  solved,  using  the  low  frequency  boundary  conditions  on  the  interfaces  to  the 
dielectric  layers.  First,  the  conductive  layer  is  divided  into  three  regions,  corresponding  to  the 
specified  boundary  data.  Using  the  separation  of  variables  techinque  on  Laplace’s  equation,  the 
potentials  within  each  region  are  determined  to  he  infinite  Fourier  series.  The  Fourier  coefficients 
are  determined  by  enforcing  boundary  conditions  at  the  interfaces  between  the  regions.  A  diagram 
of  the  conductive  layer  is  depicted  in  Figure  C.l,  where  the  x  axis  has  been  translated  by  -h,  so 
that  0  <  y  <  p,  p  =  I  —  h  For  the  resistance  calculations,  only  the  odd  symmetry  case  of 

=  — Uj  =  1/2  volt  is  considered,  since  no  electric  flux  exists  between  the  two  electrodes  with 
even  symmetry  (Ui  =  Vj). 


Figure  C.l.  Diag-am  of  resistive  layer  for  detailed  solution 


C-1 


C.l  Potential  Solution  for  Region  a 


Region  a  is  the  left  region  with  — (s  +  l)/2  <x<  s/2.  From  Figure  C.l  and  Equation  3  80, 
the  potential  within  region  a,  4>a,  must  satisfy  the  following  conditions: 


<i>a{x,Q)  =  - 


first,  we  seek  a  solution  of  the  form 


(f>a{x,  y)  =  Vaiy)  +  ^a(x,  y) 


where  Va  is  chosen  to  satisfy  the  non- homogenous  boundary  condition  at  the  lower  boundary,  A 
selected  value  of 


satisfies  Equations  C.l,  C,2,  C.3,  and  C.4.  Therefore,  Wo  must  satisfy 


Wa(x,0)  = 


(CIO) 


To  determine  we  seek  a  solution  of  the  form 


Jaix,  y)  =  'I'(x)‘^(y) 


I’sing  Equation  (Ml  in  C.7  yields 


=  0 


(C.12) 


Hence, 

<^"  xf/l' 

where  Ag  is  the  separation  constant  or  eigenvalue  for  region  a.  Thus, 
<!)"  4-  =  0 


or 


<I>(y)  =  Acos(Xay)  +  flsin(Aaj/) 
Enforcing  Equation  C.8  yields 
<I)(0)  =  A  =  0 

Thus 


(C.13) 


(C.14) 


(C.15) 


(C,16) 


(C.17) 


^(y)  =  Bsin(Aay) 

(C.18) 

<I>'(y)  =  flAaCos(Aay) 

(C.19) 

Enforcing  Eouation  C.9  yields 

(^'(p)  =  BXa  cos(Aap)  =  0 

(C.20) 

Hence 

Ag  =  —  for  n  =  1 , 3, 5, ... 

"  2p 

(C.21) 
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For  Xa  =  0,  Equation  C.15  reduces  to 

<!•"  =  0 

or 

$(y)  =  Aoy  +  Bo 
Enforcing  Equation  C.8  yields 
$(0)  =  So  =  0 
Therefore, 


(C.22) 


(C.23) 


(C.24) 


(C.25) 


=  Aoy  (C.26) 

=  ^0  (C.27) 

Enforcing  Equation  C.9  yields 

$'(0)  =  ylo  =  0  (C.28) 

which  indicates  that  Aa  =  0  is  not  a  possible  eigenvalue.  Thus,  the  complete  eigenfunction  is 

<t>r,(y)  =  Sn  sin  (C.29) 

To  determine  'i(x),  let  z  =  x  +  (s  +  l)/2,  then  '!'(?)  =  'l'(z),  and  /dx^ .  Thus 

=  0  (C.30) 


C-1 


or 


$„(?)  =  C„  COsh(Aa?)  +  «inh(Aa?) 


(C.31) 


=  CnAa  sinh(AaX)  +  D„Aa  COsh(AaX) 


(C.32) 


Enforcing  Equation  C.lU  yields  (x  =  — - 


^  implies  that  x  =  0) 


K(0)  =  E)„Aa  =  0 


(C.33) 


or  Dn  —  0.  Thus, 


'j„(x)  =  Cn  COsh(AaX) 


(C.34) 


Hence, 


=  Cn  cosh 


(2n  —  l)?r  /  s  + 


(C.35) 


Using  Equation  C.29  and  C.35  in  C.ll  yields 


CO  f* 

Wa(x,y)  -  ^e„  cosh 

n  =  l  ^ 


{2n  —  \)tv  (  s+l\]  r(2n— l)7ry 

- -  *  H - —  •  sin  - - - 

2p  V  2  ;J  L  2p 


(C.36) 


where  €„  =  B„  ■  C„.  Using  Equations  C.6  and  C.36  in  C.5  yields 

1  ^  r(2n  —  l)ir  /  s+l\l  r(2n— l)7rj/' 

y)  =  2  +  E  [-^r~  j  J  ■ — 


C.2  Potential  Solution  for  Region  b 

Region  6  is  the  middle  region  with  —s/2  <  x  <  s/2.  From  Figure  C.l,  Equation  3.79,  and 
Equation  3.80,  the  potential  within  region  b,  <jn,  must  satisfy  the  following  conditions: 


=  0 

^(^^,0)  =  0 

dy 


(C.38) 


(C.39) 
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d(j>b,  . 


0 


First,  we  seek  a  solution  of  the  form 
<Pb{x,y)  =  'I'(x)«5(y) 

Using  Equation  C.42  in  C.38  yields 

=  0 

or 

\p" 

Hence, 

_  i2 

^  ~  ~  ~  ~ 

where  Aj  is  the  separation  constant  or  eigenvalue  for  region  6.  Thus, 

=  0 

or 

<l>(j/)  =  /I  cos(A),j/)  +  Bsin(Ajt/) 

<^'(y)  =  -AAi,sin(Ajy)  +  BA),  cos{Ajy) 

Enforcing  Equation  C.39  yields 
$'(0)  =  BXb=0 

or  B  =  0.  Enforcing  Equation  C.40  yields 


(C.40) 

(C.41) 


(C.42) 


(C.43) 


(C.44) 


(C.45) 


(C.46) 


(C.47) 

(C.48) 


(C.49) 
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$'(p)  =  AXi,  sin(A4p)  =  0 


(C.50) 


Thus, 

,  nir 

Aj,  =  —  for  n  =  1, 2, 3, ... 

V 

For  A),  =  0,  Equation  C.46  reduces  to 

=  0 

or 

^(y)  =  Aoy  +  Bo 

Enforcing  Equation  C.39  or  C.40  yields 
$(0)  =  ^0  =  0 
Consequently, 

^o(p)  =  Bq 

Thus,  the  complete  eigenfunction  is 

Bo  for  n  =  0 

^(y)  = 

Bn  COS  for  n  ^  0 

Now  for  'I'(z) 

=  0 

or 

=  Cn  cosh(A6x)  +  Dn  sinh(Ajx) 
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(C.51) 


(C.52) 


(C.53) 


(C.54) 


(C.55) 


(C.56) 


(C,57) 


(C.58) 


For  Aj  =  0 


(C.59) 


(C.60) 


(C.61) 


(C.62) 


/o 

=  Bo 

■Co 

(C.63) 

90 

=  Bo 

■Do 

(C.64) 

fn 

=  B„ 

■Cn 

(C.65) 

9n 

=  Bn 

■Dn 

(C.66) 

From  Figure  C.l,  the  potential  in  the  conductive  layer  has  odd  symmetry  since  (for 
s/2  <  x  <  (s+  l)/2) 

<?i(r,0)  =  (C.67) 

<A(-x,0)  =  i  (C.68) 

Since  region  b  is  the  central  region,  its  potential  will  also  have  odd  symmetry.  Thus, 
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<i>b{x,y)  =  -4>i{-x,y) 


(C.69) 


or 

4>b(x,y)  +  4ii,(-x,y)  =  0  (C.70) 

Using  Equation  C.62  in  C.70  yields 

2go  +  2  fn  cosh  (~^)  ~  ^ 

Thus  (using  orthogonality  of  cosine  function) 

3o  =  0 

/n  =  0 

Hence, 

<^b{x,  y)  =  /oX  +  gn  sinh  cos  f  (C.74) 

n=l  \  P  /  \  P  / 

C.3  Potential  Solution  for  Region  c 

Region  c  is  the  right  region  with  s/2  <  x  <  (s  +  l)/2.  From  Figure  C.l  and  Equation  3.80, 
the  potential  within  region  c,  0c,  must  satisfy  the  following  conditions: 


= 

0 

(C.75) 

0c(r,O)  = 

1 

2 

(C.76) 

0 

(C.77) 

dx 

0 

(C.78) 

First,  we 

seek  a  solution  of  the  form 

0c(-c,  y)  =  i'c(y)  + 

■.(x,y) 

(C.79) 
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(C.72) 

(C.73) 


where  Vc  is  chosen  to  satisfy  the  non-homogenous  boundary  condition  at  the  lower  boundary.  The 
function 


Vc  —  — 
2 


satisfies  Equations  C.75,  C.76,  C.77,  and  C.78.  Therefore,  uic  must  satisfy 


du)c,  , 

if'''"’  = 


duic  /s  +  1 


ax  V  2 


For  We,  we  seek  a  solution  of  the  form 


We(x,j/)  =  ']i'(x)$(j/) 


Using  Equation  C.85  in  C.81  yields 


<|)" 


Hence, 


If  “ 


where  Ae  is  the  separation  constant  or  eigenvalue  for  region  c.  Thus, 


<I>"  +  A2<I>  =  0 


(C.80) 


(C.81) 

(C.82) 

(C.S3) 

(C.84) 


(C.85) 


(C.86) 


(C.87) 


(C.88) 


(C.89) 


O  1  O 
V-/-  itJ 


Since  this  result  identical  to  region  a,  including  boundary  conditions,  region  c  will 
identical  eigenfunction.  From  Equation  C.29 


4>„(y)  =  B„  sin 


(2n  -  l)7rj/ 


2p 


and 


(2n  —  l)7r 


for  n  =  1,2, 3, 4,... 


To  determine  ^(i),  let  x  =  x  —  (s  +  l)/2.  Then  'I'(x)  =  and  9^'I'/5x^  =  jdx^ . 


=  0 


or 


$„(x)  =  Cn  cosh(Acx)  +  Dn  sinh(Acx) 

$I,(x)  =  Cn  Aj  sinh(Aox) +  £)„  Ac  cosh(A,ric) 


Enforcing  Equation  0.84  yields  (x  =  implies  that  x  =  0) 


$;(0)  =  D„A  =  0 
or  D„  =  0.  Thus, 


il'„(x)  =  C„  cosh(A<;x) 


or 


^ri(x)  =  C„  cosh 


(2n  —  l)7r 

Tp 


Using  Equation  C.90  and  C.97  in  C.85  yields 


<^r.(x,y)  =  E  hn  cosh 


n  =  l 


(2n  -  l)7r  (  s  + 
'  X  4- 


2p 


+  l\]  .  [(2n 

-jj  ””  [- 


-  l)7ry 


2p 


have  the 

(C.90) 

(C.91) 

Thus, 

(C.92) 

(C.93) 

(C.94) 

(C.95) 

(C.96) 

(C.97) 

(C.98) 


C-ll 


where  •  C„.  Using  Equations  C.80  and  C.98  in  C,79  yields 


1  °° 

4>c{x,y)  =  + 


n=:l 


(2n  —  l)ir  f  s  +  1  y 

■{2n  -  l);ry' 

[  2p  2  )\ 

'  Sin 

2p 

(C.99) 


Referring  to  Figure  C.l,  for  s/2  <  ^  <  (s  +  l)/2 


d>(x,y)  =  4)c{x,y) 


(C.lOO) 


while 


<t>{-x,y)  =  4>a{-x,y) 


(C.lOl) 


In  addition,  for  s/2  <  x  <  (s  +  l)/2 

-^(a:.0)  =  (C.102) 

<^(-x,0)  =  i  (C.103) 


which  represents  odd  symmetry.  Thus  for  entire  problem,  we  must  also  have  odd  symmetry.  Hence, 


Hx,y)  =  -<i>{-x,y) 


(C.104) 


or 


4>c{x,y)  +  <l>a(-x,y)  =  0 


(C.105) 


Using  Equations  C.37  and  C.99  in  C.105  yields 


+  /in)  cosh 

n=  1 

or  hn  =  —e„.  Thus, 


’(2n  -  l)5r  /  s  +  1  Y 

■(2n  -  l)7ry' 

[  2p  r  2  )\ 

Sin 

2p 

=  0 


1  °° 

<l)Ax,y)  =  -  ^ 


’(2n  -  l)ir  /_  s  +  ly 

.  Cl  n 

■(2n  -  l)7ry' 

1  2p  2  )\ 

'  dill 

[  2p 

(C.106) 


(C.107) 
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C.4  Boundary  Between  Regions 


The  potential  must  be  continuous  across  the  interface  between  regions.  In  addition,  the  x 
component  of  the  electric  field  must  also  be  continuous  across  the  interface,  since  the  continuity 
has  the  same  value  in  each  region.  Due  to  symmetry,  either  interface  (between  regions  a  and  6, 
or  between  regions  b  and  c)  will  produce  the  same  boundary  equations.  Therefore,  the  interface 
between  regions  a  and  b  will  be  evaluated  (x  =  —s/2).  From  the  above  discussion,  the  region 
potentials  must  satisfy  the  following  conditions: 


d<i)a ,  S  , 
5x  ^  2  ’ 


d4>b,  S 


From  Equation  C.37 
1 

=  -  +  'Y^e„  cosh 


n  =  l 


sin 

■(2n  -  \)Ky' 

L  2p  V  2  Jl 

[  2p 

and 


d<t>a  '  (2n  —  l)7r 


dx  ^  2p 

n  =  l 

From  Equation  C.74 


e„  smh 


(2n  —  l)7r  f  s  +  1 

2^  V 


4>b  =  /o^  +  ^  Qn  sinh  cos 


f»=  I 


and 


d(t)b  ,  ^  mr  ( mrx\  ( mry\ 

=  /O  +  >  - grx  cosh  -  cos  - 

;^p  \  p  J  \  P  J 

Using  Equations  C.llO  and  C.112  in  C.108  yields 

.  oo  oo 

-  +  ^  e„r„V„(y)  =  -  ^fo-Yl  9nK^’niy) 

n=l  n=l 


sm 


(2n  -  l)7ry 


2P 


(C.108) 

(C.109) 


(C.llO) 


(C.lll) 


(C.112) 


(C.113) 


(C.IM) 
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where 


I 

■ 


I 


=  cosh 

’(2n  —  l)ff 

=  sinh 

4p 

/  MTS  \ 

K 

V  2p  y 

Kiy) 

’ 

(2n  -  l)7ry' 

=  sin 

2P 

. 

^n(y) 

=  cos  1 

'nTry\ 

K  P  ) 

(C.115) 
(CM  16) 
(C.117) 
(C.118) 


The  r  superscript  denotes  quantities  used  in  the  resistive  potential  solution.  The  /„  and  g„  terms  are 
“sifted”  by  multiplying  Equation  C.114  by  and  integrating  over  y  (0  <  y  <  p)-  Accordingly, 


n  =  1 


and 


(C.119) 


OO  CO 

n=l  n=l 

where 


(C.120) 


{/{y) 


,9iy))  =  / 

^0 


fiy)  ■  9{y)(iy 


Evaluating  the  various  inner  products  yields 


(1-1)=  r 

Jo 


dy  =  p 


and 


(^^1) 


r  ■ 

=  /  sir 

Jo 


(2n  -  l)!ry 


2p 


dy 


2p 


«-l)  = 


(2n  —  l)7r 

2p 

(2n  —  1  )7r 

2p 

(2n  —  1  )7r 


(2n  -  l)7ry 


2p 

(2n  -  l)7r 


p 

0 

-  1 


(C.121) 


(C.122) 


(C.123) 
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Evaluating  yields 


yP  f  niTy\  (mny\ 

Jo  \  P  J  \  P  J 


To  simplify  the  algebra,  let  r  =  Then  ^dr  =  dx,  and 


P 

(V'n.Vm)  -  -  cos(nr)  •  cos(mr)dr 
^  Jo 


Using  (for  integers  a  and  5) 

' 

I  (cos  az)(cos  6r)  di  =  < 
Jo 


0  for  a  ^  6 


I  for  a  =  6 


0  for  n  ^  m 


I  for  n  =  m 


Evaluating  yields 

To  simplify  the  analysis,  let  r  =  ^,  and  n  =  ~ .  Accordingly, 

’/’m)  =-  /  sin(fjr)cos(mT)dr 

Jo 


Using  (for  ^  /3^ 


(sinQx)(cos,3x)dx  =  - 


cos(a  —  0)x  cos(a  +  0)x 

2(0-0)  2(a  +  0) 


p  cos(n  —  m)r  cos(n  +  m)r]"^ 

TT  2(n  —  m)  2(n  +  m)  g 

_  p  j  cos[(n  —  m)7r]  —  1  cos[(fj  +  m)7r]  —  1 1 
V  \  2(n  -  m)  2(n  +  m)  j 

p  r  1  1  cos[(n  —  m)^]  co.s[n  +  m)7r] 


m)  2(n  +  m)  2(n  —  m)  2{n  +  m) 


Since 


cos(a  ±  /?)  =  cos  a  cos  6  T  sin  a  sin 


(C.132) 


then 


cos[(n  —  m)7r] 
cos[(n  +  m)7r] 


cos(ri7r)  cos(mT)  +  sin(nn-)  sin(m7r) 
cos(nT)  cos(m7r)  —  sin(n5r)sin(m7r) 


Using  sin(m7r)  =  0  and  cos(n7r)  =  cos  _  q  yields 


cos[(n  ±  m)^]  =  0 


Thus 


=  - - 


TT  (n^  —  rrU ) 


(C.133) 

(C.134) 


(C.135) 


(C.136) 


or 


p{2n  -  1) 


2p(2n-  1) 


2tt 


{{2r^Y  _  ^2 


Using  Equations  C.122  and  C.i23  in  C.119  yields 


TT  [(2n  —  1  )^  —  4m^] 


2p  s 

=  -7,foP 


2  •  '‘(2n  -  l)7r  2 

n=  I  '  ' 


E'-r; 


\n  =  1 


2 _ \  s  1 

(2n  -  l)7r  j  2-'°  “  2 


I'sing  Equations  C  127  and  C.137  in  C.  120  yields 


V  2p{2n-  1)  ^  p 

'‘7r[(2ri-  1)3  -4m2]]  2 


(—9m 'I'm 


(C.137) 


((M38) 


(C  139) 


{(M40) 


(^16 


(C.Ml 
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Using  Equations  C.lll  and  C.113  in  C.109  yields 
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Again,  the  superscript  r  denotes  quantities  used  in  the  resistive  potential  solution.  Sifting  Equa¬ 
tion  C.142  with  yields 
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Using  Equations  CA22  and  C.123  in  C.145  yields 
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or 
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Using  Equation  C.148  in  C.139  yields 
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Using  Equations  C,127  and  C.137  in  C.146  yields 
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Using  Equation  C.141  in  C.151  yields 
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C.5  Summary  of  Coefficient  Equations 

Enforcing  the  boundary  conditions  at  the  interface  between  regions  (for  example,  between 
regions  a  and  6)  yields  the  following  equations,  which  must  be  solved  to  determine  the  potential 
coefficients  of  the  resistive  solution: 
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where  m  is  a  positive  number  and 
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Once  the  Cn  terms  are  found,  the  terms  and  /o  can  be  found  from  the  following  relationships: 
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Appendix  D.  List  of  Symbols 


Description 

Substitution  variable  used  in  potential  solutions 

Potential  coefficients  for  layer  2 

Arbitrary  constant  used  in  potential  solutions 

Potential  coefficients  for  layer  3 

Arbitrary  constant  used  in  potential  solutions 

Velocity  of  light 

Centimeters 

Potential  coefficients  for  layer  3 

Ca^  acitance  between  electrodes 

Capacitance  from  one  electrode  to  ground 

Capacitance  associated  with  insulator 

Arbitrary  constant  used  in  potential  solutions 

Capacitance  associated  with  substrate 

Total  capacitance  between  electrodes  with  one  grounded 

Insulator  thickness,  normalized  to  the  electrode’s  width 

Potential  coefficients  for  layer  4 

Arbitrary  constant  used  in  potential  solutions 

Permittivity 

Dielectric  constant  in  layer  z  {z  =  1,2,  3, 4,  5) 

Complex  dielectric  constant  in  layer  z  (z  =  2,  3,  4,  5) 
Electric  Field 

Potential  coefficient  for  regions  a  and  c 
Farads 
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Symbol 

Kiy) 

fo 

II 

fp 

n^) 


On 


G 

9(x) 

H 

K(y} 

h 

Hz 


Im 

Pn, 


/3a 


/36 


I 

j 

j 

k 


Description 

Quantity  used  in  the  capacitive  solution 
Potential  coefficient  for  region  b 
Cricital  frequency  imposed  by  length  L 
Transition  frequency  of  the  sensor  coating 
Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Quantity  used  in  resistive  solution 
Potential  coefficient  for  region  b 
Conductance 

Quantity  used  in  capacitive  solution 
Magnetic  Field 

Quantity  used  in  capacitive  solution 

Electrode  height,  normalized  to  the  electrode’s  width 

Hertz 

Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Quantity  used  in  capacitive  solution 
Current 

Current  density 

Complex  number  equal  to  \/^ 

Wavenumber 
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Symbol  Description 


Eigenvalue  of  layer/region  z  (r  =  2,3,4,5,  a,6,c) 

L  Length  of  the  interdigitated  structure 

I  Thickness  of  sensor  coating  and  electrode,  normalized  to  the  electrode’s  width 

15  Mhos 

fi  Permeability 

n  Ohms 

w  Angular  frequency 

uiL  Critical  angular  frequency  imposed  by  lenght  L 

utp  Angular  transition  frequency  of  the  sensor  coating 

‘^3(^1  y)  Quantity  used  in  capacitive  solution 

ijja{x,y)  Quantities  used  in  resistive  solution 

Wc(x,  y)  Quantities  used  in  resistive  solution 

Quantity  used  in  potential  solutions 
Quantity  used  in  capacitive  solution 
Quantity  used  in  resistive  solution 
'I'(x)  Quantity  used  in  potential  solutions 

y^even  Quantity  used  in  charge  symmetry  analysis 

Quantity  used  in  charge  symmetry  analysis 
Quantity  used  in  charge  symmetry  analysis 
y^^odd  Quantity  used  in  charge  symmetry  analysis 

'k'  Quantity  used  in  capacitive  solution 

'kJi  Quantity  used  in  resistive  solution 

4>z{x,y)  Potential  in  layer/region  z  (z  =  2,  3, 4, 5,  a,  6,  c) 

4>^  Eigenfunction  of  layer  z  (z  =  2, 3, 4, 5) 
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Symbol  Description 


pF  Pico  (10“'^)  Farads 

p  Sensor  coating  thickness,  normalized  to  electrode’s  width 

Q  Charge  on  one  electrode 

Ql  Quantity  used  in  charge  symmetry  analysis 

Q2  Quantity  used  in  charge  symmetry  analysis 

Q3  Quantity  used  in  charge  symmetry  analysis 

Q4  Quantity  used  in  charge  symmetry  analysis 

Q5  Quantity  used  in  charge  symmetry  analysis 

Qcven  Quantity  used  in  charge  symmetry  analysis 

Quantity  used  in  charge  symmetry  analysis 

Q2  Quantity  used  in  charge  calculations 

<52  Quantity  used  in  charge  calculations 

Q2^  Quantity  used  in  charge  calculations 

<52^  Quantity  used  in  charge  calculations 

Q\  Quantity  used  in  charge  calculations 

Q3  Quantity  used  in  charge  calculations 

Quantity  used  in  charge  calculations 

Q4  Quantity  used  in  charge  calculations 

Q\  Quantity  used  in  charge  calculations 

p  Charge  density 

Rp  Resistance  through  the  sensor  coating  layer 

R,  Resistance  through  the  substrate 

cTj  Conductivity  in  layer  i  {z  —  1,2, 3, 4,  5) 

s  Electrode  separation,  normalized  to  electrode’s  width 
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Symbol  Description 


r  Substitution  variable  used  in  potential  solutions 

©'  Quantity  used  in  capacitive  solution 

Quantity  used  in  resistive  solution 
i'3(x)  Quantity  used  in  capacitive  solution 
I’a  Quantity  used  in  resistive  solution  for  region  a 

ly  Quantity  used  in  resistive  solution  for  region  c 

V'  Voltage 

Vi  Voltage  of  electrode  z  {z  =  1,2) 
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frequency  where  the  sensor  coating  changes  from  a  conductor  to  a  lossy  dielectric. 
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